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ON THE APPROACH TO EQUILIBRIUM OF SPIN 
WAVES SYSTEM. I. 


By Janusz MORKOWSKI 


Department of Ferromagnetics, Institute of Physics of the Polish Academy of Science, Poznan 
(Received May 11, 1959) 


The author calculates the spin-spin relaxation time in ferromagnetics caused by 
pseudo-dipolar interactions between spins. Calculations are based on the spin wave 
theory in the form given by Dyson. Non Hermitian Hamiltonian in the ideal model of 
spin waves system is derived and transition probabilities in first approximation are calcu- 
lated. In this approximation it appears that processes of splitting of one spin wave into 
two (and confluence of two into one) and of mutual scattering of two spin waves are the 
only possible ones. Kinetic equation taking into account these processes is derived and 
relaxation time is determined. Numerical values of relaxation time for spin waves with 
zero wave vector are calculated, showing that contribution to the line width of ferro- 
magnetic resonance absorption from such spin ‘waves is not important. 


§ 1. Introduction 

Relaxation phenomena in ferromagnetics are intimately connected with the line 
width in ferromagnetics resonance absorption experiments. Complete theoretical 
explanation of the observed line width and its temperature dependence is — in spite 
of many efforts in this direction — not yet achieved. Several relaxation mechanisms 
may contribute to the line width. The most important among them seems to be 
interactions between spin waves tending to establishment of thermodynamical equili- 
brium in spins system (Akhiezer 1946; Kasuya 1954 a; Kaganov, Tsukernik 1958), 
interactions between “ferromagnetic“ d-electrons and conduction electrons (Abra- 
hams 1955, Mitchell 1957, Kasuya 1954 b, Kondoh 1953), and interactions of spin 
waves with lattice vibrations (“phonons”) levelling spin system and lattice tempera- 
tures (Akhiezer 1946; Kittel Abrahams 1953; Sugihara 1953, Terasaki, Manari 1957; 
Kaganov, Tsukernik 1959 also Akhiezer et al. 1959). Another attempts ie explain line 
width were also published (see for instance recent concise review or this question 


The situation is thus quite complicated, relative importance of these various 
s is not exactly known. It seems thus that 


(5) 
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possibilities of improving the existing theoretical values of relaxauon times by mathe- 
matically more rigorous treatment cannot be overlooked. 

The purpose of the investigations presented here is to reexamine the theory of 
relaxation processes in spin waves system. Existing papers on the same problem 
(Akhiezer 1946, Kasuya 1954, Kaganov, Tsukernik 1958) are all based on Holstein- 
Primakoff’s theory of spin wave interactions. Recently Dyson (1956 a, b) critized the 
formalism of Holstein and Primakoff shoving that: 1) spin-wave amplitude operators 
introduced by them do not have proper physical interpretation and 2) mathematical 
approximations indispensable for simplification of intricate Holstein — Primakoff's 
Hamiltonian are not always justified. In Dyson’s papers more exact theory of spin 
wave interactions is given. In the present paper ‘we take Dyson’s treatment of spin 
waves as a starting point for calculation of relaxation times, hoping that values ob- 
tained in this way are more appropriate than the former ones. 


§ 2. Hamiltonian 


We shall discuss the relaxation of spin waves in perfect cubic crystals of a fer- 
romagnetics. We start with the Heisenberg model of ferromagnetics. The s—d inter- 
actions will be completely neglected, although it is well known that the contribution 
to the line width caused by conduction electrons is by no means negligible in real 
ferromagnetic metals. Our task is merely to examine relaxation processes connected 
with direct interactions of spin waves. 

We assume that each lattice point is determined by a lattice vector ۳ let N 
be the total number of lattice sites in the crystal. Periodic boundary conditions are 
assumed. At each lattice point j is situated spin vector رگ‎ = (Sf, S?, S?) of magni- 
tude S. 

The Hamiltonian of the system consists of 3 terms: 


1) Zeeman energy in a uniform external magnetic field H, directed along the 
z-axis of a frame of reference: 


~ 


2ug H DS | (2.1) 
3 


Hg — Bohr magneton; — 
2) isotropic exchange interaction with (positive) coupling constant Jin (in practice ` 
AJ = J, when j j and k are nearest neighbour lattice sites and J, = 0 for other cases) 


- و‎ NSS (2.2) 


3) anisotropic interaction energy between spins. We adopt the lowest order aniso- 
tropic coupling i. e. the so called “pseudo- dipolar Fees of the form (see: Vans : 
Vleck 1937, Van Kranendonk, Van Vleck, 1958): à 


X Da [S 8, - 3 ° o ez? (2.3) 


=. 


N 


On the Approach to Equilibrium of Spin Waves System 5 


wherein 
eg = (j HU — H 
The pseudo-dipolar coupling arises from indirect effect of spin-orbit interaction; the 
coupling constant D,, rapidly decreases with increasing distance between j and É, 
being significant only for nearest neighbour spins. The classical magnetic spin-—spin 
interactions are also of the form (2.3), the coupling constant for these interactions is 
Dj, = (2ug*/|j — E (2.4) 
They are weaker than the former ones at nearest neighbours distance but have a long 
range character. For S — 1/2 higher multipole interactions can also be important, 


but we shall neglect them in this paper. 
The total Hamiltonian is thus equal to 


J J+k j+k 
where 
L-—2ugH (2.6) 
In the following considerations we shall often use the well known operators 
SF = S¥ + iS? 
Ds ne o Q3) 
and their Fourier transforms: ۱ 


EN SC SE 
J 


eg (2.8)‏ و 
j‏ 


§ 3. Formulation of the problem 


We shall examine the processes of mutual scattering of spin waves, leading to 
establishment of a thermodynamical equilibrium state of spin system and we shall 


find the spin-spin relaxation time. y 
As the basis of our calculations we adopt Dysons (1956 a, b) theory of spin wave 


interactions. E i 
According to Dyson (1956 a) we define the spin wave state containing ره‎ spin 


waves with wave vector À by: 
la > = H (2S) ^4 (a5) * (S2) 2| 0 > (3.1) 
a 


where the ground state fulfils the condition: 


d .ود‎ | 0 < =0 (3.2) 
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The states | a > (with the exception of case Y a, = 1) are not orthogonal to each 
n 


other. The spin waves are thus not independent, some constraints are imposed on 
them — the so called by Dyson kinematical interactions. 

In order to avoid mathematical difficulties connected with kinematical interac- 
tions and manifesting themselves in nonorthogonality of states | a > Dyson intro- 
duces another complete set of orthonormal states | a), being in one-to-one correspon- 
dence with the primary set of states | a >. The “ideal spin wave states“ | a) are 
defined by: 


la) = IT (a,!)~* (aj )**| 0) (3.3) 
A 


2,10) = 0 (3.4) 


where aj, ره‎ are the well known harmonic oscillator operators (creation aj and anihi- 
lation a, operators), satisfying the commutation rules: 


EI aj] = Oru [æ;, a] = la}, ai] = 0 (3.5) 
and 
a; a,|a) = azla) (3.6) 
or formally 
ai = a, 


Starting with the Hamiltonian # (2.5), expressed in terms of operators و رگ‎ Sy, SF 


we can — following the method given by Dyson — construct an operator # in the 
ideal spin waves space (“ideal model“), operating on the ideal states | a) in the same 
way as does Hamiltonian # on the physical states [a 

The analysis given below leads to an other simple formal way of constructing 7 
more useful than Dyson’s original one in cases of complicated Hamiltonian. 

Let 2 be an Hermitian operator in the physical model and | y > any physical 
state. In the ideal model the state | y > is replaced by the ideal state | v). The state 
(v | adjoined to | y) in the ideal model is defined according Dyson (1956 b) by (y | — 
(y |F where the operator F is the so called "indefinite metric tensor". The cor- 
respondence between matrix elements of an operator Q in the physical and ideal 
models is given by equation (Dyson, 1956 b): 


> 9|9|v > = (ell y) = (p |F y) | (3.7) 


(symbol ~ means, that the quantity in question refers to ideal model). It is easy to 
verify, that operator 2, Hermitian in the physical model, satisfies in the ideal model 
the condition (Dyson, 1956 b) 


FQ = QF (3.8) 
(œ denotes the Hermitian adjoint of the operator 2) i. e. Q is self-adjoint with 


respect to metric F. 
: i 
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The explicit expression for the metric tensor F can be given using a particular 


set of states. Let us take a set of integers 0 < u; X; 2S, corresponding to lattice points j. 


We define, according to Dyson (1956 a, eq. (30), a complete and orthogonal (but not 
normalized) set of states with localized spin deviations: 


lu > = II (28) +" (uj) 5 (S+)4| 0 > (3.9) 


In the ideal model the state | u > is connected with the ideal state 
hu) = IT (uj) * bei 0 (3.10) 
j 


describing a collection of harmonic oscillators attached to lattice points 7, where 
oscillator amplitudes 7° (creation operators) and 7; (annihilation operators) satisfy 
commutation rules 


fo, gel än m md = m.m] = 0 (3.11) 


and the integer u, is the eigenvalue of operator 5j n;. 
By (3.7) we have (Dyson, 1956 b): 


(wiF |v) = < ulv > (3.12) 


From the other hand from (3.9) follows (Dyson, 1956 a, eq. (31)) 


le (3213)‏ ری < < و | هه > 

where à, = II SCH (du, e is the Kronecker symbol) and 
` 
: 1 ولا‎ — 1 1 

From eqs. (3.12), (3.13) it follows that: 

(ul F |v) = 
This is equivalent to: ۱ 

(u| F = F, (ul (3.15) 


It is not hard to prove now, that | 
1 "5 à 
GE, = (ol UT vj? (TT ۲ = eu | (1- Lun) n| Dit, (3.16) 
j 


It is thus evident, that correspondence 
lu > > lu) 
j <> = (o| F (3.17) 


4 
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can be realized by replacing operators S}, S7, S by n*,n; according to the rules 
(c. f. Maleev, 1957): 


= N 
رگ‎ > GSM 0 SEC nj 1 N 
S; > خ(ک2)‎ nt (3.18) 
SF = n Ces 5 
Replacing operators $7,.57, S in Hamiltonian # (2.1) by n,n; we formally 


get the operator Z in the ideal model. To pass over from localized oscillators to spin 
waves the following Fourier transformations must be performed: 


nj = Nê 2 e qt 
ny NI D elg, (3.19) 


It is easy to prove, that (3.18) and (3.19) are equivalent to 
رد‎ (25) af 


$7, (25) [a, CSN SI oi, ;a, ay) (3.20) 


In transforming Hamiltonian (2.1) by means of (3.20) some approximations are 
used which are discussed below. Sums of the form 


N^! E (äs e P E X) 
J+k 
may be written as 
NO y e id E Dy e-i G- 5 

k IGEk) 
The pseudo-dipolar coupling constant is a function of h = | SÉ H only, rapidly 
decreasing with an increase of h. If we neglect the lattice points lying on the ends of 
the crystal the sum X D, exp [—iA- (j — k)] does not depend! on E (compare Hol- 

J 
stein, Primakoff, 1940; page 1102). 
Thus we have approximately: 


News D eO AE) A (A Së u) 22 D(h) eg ® (3.21) 
ve ییحی‎ E CD MM . cic 


1 The formal grounds for this approximation lies in assumption of periodic boundary conditions. 
On the other hand theoretical and experimental evidence points out that such approximations must be 
used with some care — compare for instance Kittel 1958. 
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After treating in this way other sums and after rearranging the elements the 
quasi-Hamiltonian takes the form: 


y — z RE = "Ds ł 
H = N (PGS IS) RE E Dat mt D Qf af at 


Lu ee A @SNM ر2‎ (Tr aj at ya, + Tt aa, 82) — 


20 
— (25N)-1 3 ) + Ry SP Jah. a a,a, — (3.22) 
مب‎ p Sort =; 
(SN) 2 Qj وبوبره‎ 4, €, Lo + (2SN*) À 2 TT رت‎ ra Pe 000 
E (25N)-? 2 Qi qi «€, a, Ray Zo c, 
QO 


where the coefficients are defined by 
Qi = 4S (Pi — Py + 2iT?) 
RE = 4S (Pj + Py + 2TP) (3.23) 
FE Te wh 
wherein 


Pê = - 3 D [J (h) — D (h) + 3D (AA? oF 


LE 
= &mn-x,y,2 


TP — 3 Y D (h) (hg hh) e À * ` (824) 
hz 0 


Summations over h extend on all neighbours of the lattice point j = 0, practically 
we can, however, consider only the contributions coming from nearest neighbours 


and neglect others. 
It is evident from (3.22) that the quasi-Hamiltonian is not Hermitian. The method 


of calculation using this unusual quantity will be explained below. 
The quasi-Hamiltonian # may be divided into *nearly" diagonal part #, and 


interaction Hamiltonian €, which will be of primary importance in our calculations. 
(RES NITET TU ^ (3.25) 

where E, = N (PS? — LS) and 
X = F(R} + Ry — 2SP§ + 1) aj a, + D (Qi af ats + Qi me) (826) 

n n 

The coefficients R} + Ry — 2SP6 are in first approximation proportional to the 
exchange integral between nearest neighbours J, while Q, are proportional to the 
pseudo-dipolar coupling constant D. As D|J <1 the Q; term ıs a small correction. 


We are interested in energy of spin waves only in first approximation thus the correc- 
tion terms will be neglected (for more detailed discussion see Appendix I). 
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In this approximation we have the energy €, of spin wave with wave vector A: 


Ri + Ry - 258 هراد‎ 5 9 J()(lL et) i ae, (3.27) 
hZ0 - 


In nearest neighbours approximation we get 


e, = L4 JS (yo — Ya) (3.28) 
wherein 
n= 3 (3.29) 
(a 


and y, is the coordination number of the lattice. 


8 4. Transition probabilities in the ideal model in first approximation 
Let |y > be a time dependent state in the physical model, satisfying the Schrö- 
dinger equation 
9 
ih — |y > = H| y > (4.1) 
ot 
The corresponding equation in the ideal model is (compare Dyson 1956 b, 82) 


2 ae 
ih = ly) = Hl y) (4.2) 


where |y) is an ideal state corresponding to physical state |y >. The equation 
adjoined with respect to the metric F to eq. (4.2) is 


P P 
ihn (yl F = ) ۶ (4.3) 


where we have taken into account the eq. (3.8). i 

Equation (4.2) determining the vector of state |y) and eq. (4.3) determining the 
adjoined vector of state (y| F form the starting point for our calculation of transition 
probability from one ideal state |a) (defined by eq. (3.3)) to another one. 

The derivation of transition probability is well known and we quote here some 
important points of this calculations merely because some modifications are needed 
here in order to account for indefinite metric and non-Hermitian Hamiltonian. 

We expand the state vectors |), (y| F into orthonormal set of states la): 


ly) = X la) (aly) 


el P = Y (vl Fla) (a (44) 
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with expansion coefficients depending on time. The probability that at some given 
time the system is in the state |a) is defined as: 


(pl F |a) (a |y) (4.5) 
From eqs (4.2), (4.3) it follows that the quantity (y| F ly) = X (y| F |a) (a |y) has 


a stationary value and thus can be normalized to unity. 


According to eq. (3.25) the quasi-Hamiltonian # (apart from constant Ej) may 
be broken into diagonal (and Hermitian) part X, and interactions part x y Denoting 
by e, the unperturbed energy of the state |a), given by the Schrödinger equation 
(compare (4.2)) with disregard of the €, 


Eq |a) = Lola) (4.6) 
and introducing other expansion coefficients G (t), C, (t) by 


(aly) = C, () e & *' 


(yl Fla) ent (4.7) 
we obtain from eqs (4.4), (4.2) and (4.3) 
SCH Lxx z 2 (a | Sri b) Cs (t) ef TIU 
> 


(4.8) 


dC, (t j SN ya = 
۵ - ۶ ۶۳ ۰ ۳606۱3۸9 
b 
As is well known equations (4.8) can be solved in first order approximation of the 


perturbation theory in the usual way (see for instance Schiff 1955). As a result we get 
the transition probability per unit of time from the ideal state |a) to the state |b) 


Wa = 7 (a ار‎ 8) bl dlame) (49) 


It is worth noting here that رم‎ evidently secures fulfilment of the principle 
of detailed balance in spite of non-Hermitian Hamiltonian #. In this case detailed 
balance follows from eq. (3.8), which replaces Hermicity condition of the traditional 


derivation. 


8 5. Kinetic equation and relaxation time 


We proceed now to the analysis of matrix elements involved in the expression 
(4.9) for transition probability. In the first place it is not difficult to see that some 
contributions to the transition probability IP vanish identifically for all states 
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|a) and |b) 
(al D Qy alor ده‎ €, Me lb) (61 DB) OF rere %4 Zo La la) = 0 (5.1a) 
Ago Aga 
(ol D Ti ak, 44 E I) (bia aa, at ,@,|a)=0 (5.1b) 
Aug Aug 
(al >) Qi CSC +A Z, a, CES Go Q, |b) (b| 22 Qi a Gy, a, 2 aE 
go voa 
X a, Xy |a) = 0 (5.1c) 


Thus in the approximation adopted here establishment of thermodynamic equilibrium 
in spin-wave system is secured by processes involving only 3 and 4 spin waves. There- 
fore, the interaction part € ر‎ of the ideal Hamiltonian, determined by eqs (3.22) — 
(3.26) is in the present calculation equivalent to 


X =H + X, (5.2). 
where 
Ha = CSN (Py ap OC Ta) (5.3) 
20 
H, = — QSN)3 Y (RF ER, OSP at aies (5.4) 
Aoc 


a) Transition probabilities for processes with 3 spin waves 


It is convenient to transform #, of eq. (5.3) into more symmetrical form 
E — y+ gt + + 
H = 2 (d10 %1 Ge Lite + Aie ره‎ % Xito) (5.5) 


where 


A$ = — } (2SN WI (T£ + T#) (5.6) 
Transition probabilities leading to creation of spin wave with wave vector A are: 
(aur ap مه‎ > (apt ay ag = On] A) A AZ, (a, + 1) a, a, Ô (e, — €, — 
| - e-7 -7( 
¥ A ap gel + (a+ 1+ 1, ag—1) 7 (8x/h) A5, A, (a, + 1) (a, + 1) a, 9 (e, + 
+e,—e)ô(i + —9) ^ (5.7) 
W 3 45,45) (a1 a-1,2,11) = (87/ h) A, Aj, (a, + 1) a, (a, + 1) 6 (e, — e, + 
+ 6) (à — » + 0) 
where €, stands for the spin wave energy determined by (3.27) and 
80 21 for À 20 
= 0 for ÂZ0 (5.8) ` 
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Transition probabilities for annihilation of spin wave u can be easily obtained from 
(5.7) by suitable interchange of indices, for instance: 


(3) = 
Wa, dy; ag) + (@4—1, a, 1, 454-1) = (8r/ h) AZ Zi a, (a, + 1) (a, 1) ö (€, at; 


— £) ô (u — ? — o) (5.9) 
b) Transition probabilities for processes with 4 spin waves 


Ha may be written as 


Hy = 2 CT a aa (5.10) 
where (5.11) 
— 4N Lë 3 J (y, ch نج‎ EE Vota) == (2A, T 2A, وم‎ ar A, ge A+) 
A, =} J (l — 3A2/h?) D (h) exp (— À - h) (5.12) 
h+0 


. Transition probability per unit of time for creation of spin wave li is 
we, a,» ay, G5) (au tl, ay را‎ ay- ag—1) ~ (Ba hy Va Y oar) CY E LDS 
Ra, repas El) a 59 (6 GG, GC e)à(u4-x—»— 6) (5.13) 
۱ 
The other two expressions which are to be considered can be found by cyclic inter- 
change of indices. 
For annihilation on the spin wave D we have 


(4) ae) dot 


(Ay ay ay, a,)+(a,—1, a4—1, ay+1, ag = (8x) h) ( 
Sé a, a, (a, + 1) (a, + 1) 9 (e, F Ex — ره‎ — €o) Ô (à + x—»— 9) (5.14) 


and other two expressions by interchange of indices. 

In order to get the kinetic equation from transition probabilities given above 
average numbers of transitions must be obtained. The process of averaging transition 
probabilities over suitable ensemble of systems at given conditions, simple and rs- 
quiring no comments in the traditional approach, must be discussed in some detail 
in our case since in this averaging kinematical interactions must be taken into account. 
We can simplify the discussion by assumption (implicitely involved in already cited 
papers of Akhieser 1946, Kaganov, Tsukernik 1958) that average of the product of 
quantities a, is approximately equal to the product of averages. This assumption 
seems to be quite satisfactorily fulfilled because the dynamical interactions between 
spin waves are weak. The problem is thus reduced to obtaining the average values of a, 
at actual state of the system taking into consideration. kinematical interactions: 

When the spin wave system is not far from thermodynamical equilibrium it may 
` be assumed (compare for instance Akhieser 1946, Kaganov, Tsukernik 1958) that 
the average occupation numbers are equal 


, à,— وق‎ 5 x 
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where à? is the average equilibrium value of a, at given temperature of the spin-system 
and the time dependent part £, is small compared with aj. In most experimental 
situations this assumption seems to be well fulfilled excluding only ferromagnetic 
resonance saturation experiments at high microwave power. 

Since, as it shall be evident from calculations given below, the quantities é, do not 
enter into the formula for relaxation time, only average values at equilibrium aj 
must be determined. Fortunately, this intricate problem of evaluating à? when kine- 
matical interactions are taken into account is in principle solved in Dyson's paper 
(1956 b). Dyson shows namely, that at temperatures low compared with the Curie 
temperature of ferromagnetic in question, say for 0 < T < 1 T, the contributions 
from non-physical parts of ideal states to the partition function are negligible. In other 
words it may be said that in evaluating the partition function of the spin-waves 
system the restrictions imposed by existence of kinematical interactions are not import- 
ant. We can thus neglect them in practice. Thus we get for average equilibrium values 
of occupation numbers à (neglecting in this point weak dynamical interactions bet- 
ween spin waves) the well known formula: 


ag = (ef — 1y3 (5.16) 
B = (kT) 


The mean net rate of change of the number a, of spin waves with wave vector A 
is calculated as the mean difference between the number of spin waves created as 
a result of collisions per unit of time and the number of annihilated ones per unit of 
time. The resulting kinetic equation is: 


(a) QUU Gn I AZ AL UG, + 1) a, à, — à, (2, + 1) (& + 1] ۵ (e, — e, 


— e) ô (ù — 7 — 0) + 24% Am [(à, + 1) (à, + 1) à, — à, à, (à, + 1)] 8 (e, + £, — 


— &) ۵ (4 + » + (Arf) DIE + VUN HIT + 1) (2, + 
+ 1( 4,4, — à,à, (à, + 1) (à, + 1)] ê (e, + e, —e,— e) ô (H+ — v— o) (5.17 


where we have taken into account the symmetry properties of coefficients je du 
Putting رت‎ given by (5. 15) into the equation (5.17) and — on account of |£ N <a; — 
retaining only terms linear in deviations €, the eq. (5.17) can be reduced to the form 


day = 1 m 
ak =- be 3 S . (5.18) 


In deriving this equation the formula (5.16) was taken. into account and that thus 
expressions of the form , 


]) + 1) a; a — ā (a2 + 1) e 1 mm ô e m8 سح‎ SES 0 | (5.19) 
vanish. ` | 
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Equation (5.18) defines relaxation time v7 for spin waves ji, (r+)-lis thus given by: 
(a) = Gr) D (47, An (2; + 80 D 8 (e, — &, — (ه‎ 8 — 5 — 8) + 
+ 24,4, — à) ۵ (e, + e, — e) 6 (u + Y» — 6)} + 
+ Gel h) IE HEINE + و۵‎ + 1) — 


— GC a9] ô (e, +e, — e, — e) Ô (un + x — v — 6) (5.20) 


This equation forms the basis for calculation of relaxation time. A case important 
for practical purposes is discussed in $ 6. 


8 6. Relaxation time for spin waves with zero wave vector 


As a simple application of the above consideration we will discuss practically 
important case i. e. relaxation of spin waves with wave vector u = 0, which are first 
of all excited in conditions of microwave resonance absorption experiments. It is 
expected (Kasuya 1954, V. Kranendonk, V. Vleck 1958) that as a result of several 
mechanisms spin waves with nonzero wave vectors may be also excited in ferromagne- 
tic resonance. More detailed analysis of this question will be given in subsequent paper. 

As usual the analysis will be limited to long spin waves only, i. e. we assume 
dû << 1 where 6 is a vector between neighbouring lattice points. In this approxima- 
tion the energy of spin waves takes the well known quadratic character: 


e, 2 L + oi? 
w = JSa? (6.1) 


(a is the lattice constant). 

It is evident, that for spin wave = 0 transition probabilities (5.7) for 3-wave 
processes are vanishing, as the momentum and energy conservation laws (expressed 
by 6 — factors) cannot be fulfilled. Only 4-wave processes have to be discussed. 

For pseudo-dipolar interactions coefficient A, (5.12) may be evaluated in 
nearest neighbour approximation with coupling constant D. The result for long spin 
waves is 


A, — p (4? — 34) (6.2) 
with , 
p. = $ Da? for simple cubic lattice 
Pree = 0 for body centered cubic lattice . (63) 


Pree = À Da? for face centered cubic lattice 
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Inverse relaxation time is thus obtained from eq. (5.20) for waves with # = 0 | 


TL (270 p*] hoy Ne 302 + 0?) — 2», ol 


x [v2 + شم‎ — 3 (92 + o3] [22 (87 + ag Fl) — & à; (6.4) 
x ê (o) ô (x — v— o) 
The sums can be evaluated by transforming them into integrals in the usual way: 


316...) = N (aln) 27)? f d» (...) (6.5) 


where the numerical factor n equals 1, 2, 4 for simple cubic, body centered cubic 
and face centered cubic lattice respectively. As the main contributions to sums in 
(6.4) come from small values of wave vectors the integrals may be extended to the whole 
reciprocal lattice rather than confined to the first Brillouin zone only. Introducing 
polar coordinates, integrating over polar angles and passing to new integration va- 
riables: 

x = Boc? 


y — Bo»? 
we get finally: 


E 33 DNA EER BE 
v "agn (7) iss) | af + 
0 0 


1 1 1 1 1 
x EE ۳ ۳ T Bem 
wherein: 


e = BL = 2ug H]kT ۱ (6.7) 


and r equals 1, 0, 2* for the simple, body centered, and face centered lattices, respecti- 
vely. 

In the case e <1 (fer weak magnetic fields and not very low temperatures) the 
inverse relaxation time varies approximately as T4. For very low temperatures £2» 1 
and و۲‎ is proportional to T* exp (—2ug H/kT), rapidly tending to zero as the tempe- 
rature tends to zero. 

As an example the numerical values for nickel may be given. 
For Ni: 
Tfec = 24 
| vor 
J 22-230k-1°K (see Kittel, 1954); the factor 2 appears owing to different defini- 
tions of exchange integral in Kittel’s book and in the present paper. For £ <1 we get 


79! 2 1.8 x 109+ (D/J)? (T/230°)4 | (6.8) — 


E. 
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The ratio of pseudo-dipolar coupling constant D to the exchange integral J is very 


difficult to estimate. As a first, very rough estimation it may be taken as (Van Vleck 


1956) 
DIJ ~ (g -2)? 


where g is an experimental value of the spectroscopic splitting factor. From ferro- 
magnetic resonance measurements it follows (Kittel 1954) 


g 2 2 
and so from (6.8) 
To! ~ 108 (7/230°)4 sec-1 (6.9) 


The order of magnitude of the relaxation time indicates that the processes consi- 
dered here cannot determine the line width of ferromagnetic resonance. 

Comparing our result with the results obtained by Kasuya (1954 a) on the basis of 
Holstein —Primakoff's theory we may remark that: 1) in our considerations exchange 
forces do not contribute to relaxation of spin waves with zero wave vector, 2) the 
relaxation time varies in our case as T-* (for not very low temperature) while Kasuya 
found 7-3. For nickel the value of 7, at room temperature obtained from (6.9) is 
several times shorter than value which may be deduced from formula (3) of Kasuya’s 
paper (1954 a). 

To complete the considerations presented here relaxation time of spin waves 
with nonzero wave vectors must be calculated. This will be done in the second part 
of this paper. 

I wish to thank Professor S. Szczeniowski for encouragement and valuable dis- 


cussions. 


Appendix 1. 
It may be shown that omitting of the terms 


22 (Qi ar at; +0, eg À (A I.1) 


a 


in Ho of eq. (3.26) is equivalent to the neglect of quantities of the order of (D/J)? 
as compared to unity (in connection with this section compare Akhiezer 1946). By 
means of Holstein — Primakoff’s (1940, eqs. (16)—(18)) transformation 


guo beni ib bag 


a_, = 2-3 ei-1(b, — b_,) (A L2) 
i b =; cr ly SÉ 
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where g, is polar angle of the vector A, 
d ea + (eå (۳ 
"2 ee 2O 
and (A 1.3) 


p SEET 
dE 


IO)‏ = ارم 
of eq. (3.26) may be diagonalized‏ 7 
Hy = De — 20,9% ef e, (A L4)‏ 
1 


and interaction Hamiltonian #, expressed in terms of operators c7, c}. In the nearest 
neighbours approximation and for long waves the ratio: 


2Q,l/e, ~ DIJ (A 1.5) 
is small as compared to I. Thus approximately, with error of order (D/J)? 
| (2 — DO c e, (A L6) 
121 1420 
and (A I.4) reduces to: 
#2 2 acte (A 1.4’) 


Returning to the old operators aj, a, eq. (ALA) is replaced by: 
Ho 2 2( ره‎ af a, (A 1.7) 
a 


and interaction Hamiltonian takes the primary form determined by equs (3.22)— 
(3.26). It is thus proved that. omission of operator (A I.1) is equivalent to approxima- 
tion (A I. 6), connected with rather small error as the condition (D/J)* <1 is-in 
practice well fulfilled. 
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ÜBER EINE EXPERIMENTELLE METHODE ZUR 
UNTERSCHEIDUNG DER VERSCHIEDENEN ARTEN DER 
MULTIPOLSTRAHLUNG IM RÖNTGENSPEKTRUM 


I. VERALLGEMEINERUNG DES REZIPROZITÄTSTHEOREMS FÜR 
ELEKTROMAGNETISCHE MULTIPOLSTRAHLUNGSQUELLEN 


von A. RUBINOWICZ 


Physikalisches Institut der Polnischen Akademie der Wissenschaften 
Institut für theoretische Physik der Warschauer Universität 
(Eingegangen am 15 Mai, 1959) 


Das von H. A. Lorentz im Jahre 1896 für die elektrische und magnetische Dipol- 
strahlung angegebene elektromagnetische Reziprozitätstheorem wird für elektrische und 
magnetische Multipolstrahlungen beliebiger Ordnung verallgemeinert. Für den Fall der 
elektrischen und magnetischen Dipul-sowie der elektrischen Quadrupolstrahlung 
werden Formeln angegeben, die die im verallgemeinerten Reziprozitätstheorem auf- 
tretenden Entwicklungskoeffizienten durch die elektromagnetischen Feldstärken und 
ihre Ableitungen ausdrücken. Das Theorem ist sowohl zur Lösung technischer Probleme 
(Ausbreitung von Radiowellen von Antennen endlicher Größe und beliebiger Gestalt) 
als auch physikalischer Probleme (Bestimmung der Art der Multipolquellen im Gebiet 
der Röntgenstrahlung; Lösung exakter elektromagnetischer Beugungsprobleme für 
beliebige Multipolquellen, wenn solche für elektrische und magnetische Dipolquellen 
bekannt sind) verwendbar. 2 


Die Frage, welcher elektrischen oder magnetischen Multipolstrahlung eine 
Spektrallinie ihre Entstehung verdankt, kann im optischen Wellenlängenbereich 
des Zeemaneffektes entschieden werden. Im Gebiete der 
Röntgenstrahlung versagt dieses Hilfsmittel jedoch vollkommen. Und gerade hier 
wäre es interessant eine Methode zur Ermittlung des Multipolcharakters der Spek- 
trallinien zu besitzen, da ja im Röntgenspektrum mit wachsender Atomzahl die Strah- 
lung der höheren Multipole gegenüber der elektrischen Dipolstrahlung immer mehr 
hervortritt. Bei Uran beträgt z. B. die Intensität der elektrischen Quadrupollinien 


bereits ungefähr 1/10 der stärksten elektrischen Dipollinien (Rubinowiez-Blaton 1932). 


Eine andere Methode zur experimentellen Ermittlung des Multipolcharakters 
Weitwinkelinterferenzen. Bereits im 


der elektromagnetischen Strahlung bieten die 
(21) 


experimentell mit Hilfe 
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Jahre 1911 hat P., Selenyi (1911) gezeigt, daß sie zum Nachweis der elektrischen 
Dipolstrahlung im optischen Wellenlängengebiet verwendet werden können. Aber 
auch in dem uns interessierenden Róntgenstrahlengebiet treten sie in Erscheinung, 
nümlich bei den Kosselschen Versuchen (Kossel, Loeck und Voges 1935, Kossel 
1937) mit Róntgeninterferenzen an Gitterquellen in Kristallen. Ich habe daher schon 
vor einigen Jahren (Rubinowicz 1949) darauf aufmerksam gemacht, daD diese Inter- 
ferenzerscheinungen zur experimentellen Ermittlung des Multipolcharakters der 
Röntgenstrahlung der Atome benützt werden könnten. Dazu wäre nur die experi- 
mentelle Bestimmung der relativen Intensitäten der Kosselschen Beugungsstreifen 
erforderlich. Allerdings ist die in Rede stehende Erscheinung sehr lichtschwach und 
wurde auf photogräphischem Wege bisher nur bei den stärksten elektrischen Dipol- 
linien (K, und K,) beobachtet. Da jedoch die begründete Aussicht besteht, daß zur 
Beobachtung der Kosselschen „vollständigen Reflexsysteme“ die ungleich empfindli- 
cheren und für Intensitätsmessungen sehr bequemen Photoelektronenvervielfacher 
verwendet werden können, glaube ich daß es bereits an der Zeit ist, eine entsprechende 
Theorie für die elektrischen Quadrupol-Röntgenlinien zu entwickeln. Dies ist umso 
notwendiger, als diese Erscheinung wenigstens gegenwärtig unsere einzige Hoffnung 
ist, auf experimentellem Wege den Multipolcharakter der Röntgenstrahlung der 
Atome und eventuell auch der hinreichend weichen y-Strahlung radioaktiver Kerne 
zu bestimmen. Allerdings muß bemerkt werden, daß selbst bei den elektrischen 
Dipollinien einige wenige experimentelle Einzelheiten noch nicht vollständig geklärt 
sind (vgl. Kossel 1937, v. Laue 1937, 1941, 1948). M. v. Laue (1937, 1941, 1948) ist je- 
doch der Ansicht, daß eine bessere Näherung seiner dynamischen Theorie der Rönt- 
geninterferenzen zu einem Verstündnis der bisher ungeklürten Einzelheiten führen 
muf. 

In der vorliegenden Abhandlung soll eine Vorarbeit zur Inangriffnahme des oben 
formulierten Problems erledigt werden. Max von Laue (1935, 1937, 1941, 1948) 
hat bei der elektrischen Dipolstrahlung einen sehr eleganten Weg eingeschlagen, 
der auf der Verwendung des Reziprozitütstheorems! für die elektrische Dipolstrah- 
lung beruht. Auch in dem Falle der elektrischen Quadrupolstrahlung ist dieser Weg 
gangbar, wenn nur das Reziprozitätstheorem entsprechend verallgemeinert wird. 
Eine solche Verallgemeinerung und zwar sogleich für die elektriséhen und magneti- 
schen Multipolstrahlungen beliebiger Ordnungen sollen die nachstehenden Betrach- 
tungen liefern. Sie scheint mir nämlich auch für andere physikalische und technische 
Probleme, bei denen Multipolstrahlungen auftreten, von Interesse zu sein. 

In der Physik kann man z. B. aus den strengen Lósungen der Beugungsprobleme 
für beliebig im Raume orientierte elektrische und magnetische Dipole exakte Lósungen 


der gleichen Probleme für beliebige elektrische und magnetische Multipolquellen 
herstellen. | 


! Der erste Beweis für das elektromagnetische Reziprozitätstheorem für die elektrische und magne- 


tische Dipolstrahlung rührt wohl von H. A. Lorentz (1896) her (vgl. dazu Rubinowicz 1957, S. 242). 
Für diese beiden Fälle hat auch Sommerfeld (1931) einen Beweis angegeben. z 


en w 
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Als ein technisches Problem sei das Problem der Ausbreitung der Radio- 
wellen genannt, die von einer Antenne endlicher Größe und beliebiger Gestalt emit- 
tiert werden. Die Strahlung einer solchen Antenne kann man nämlich durch eine 
geeignete Superposition von elektrischen und magnetischen Multipolstrahlungen 
erhalten. Da man mit Hilfe des Reziprozitätstheorems elektromagnetische Felder 
höherer Multipolquellen aus denen einer elektrischen und einer magnetischen Dipol- 
quelle erhalten kann, ist somit grundsätzlich die Lösung des Problems der Strahlung 


‘einer beliebigen Antenne aus der Lösung des gleichen Problems für Dipolantennen 


ableitbar und zwar sowohl im Falle der ebenen als auch der kugelförmigen Erde. 

Falls es künftig gelingen sollte, die Kosselschen Röntgeninterferenzen zu beo- 
bachten, die von der y-Strahlung radioaktiver Kerne erzeugt werden, könnte man 
das in der vorliegenden Arbeit angegebene Reziprozitätstheorem zur Berechnung 
der für die einzelnen Multipolstrahlungen höherer Ordnung charakteristischen Inten- 
sitäten der Kosselschen Interferenzstreifen verwenden. 


8 1. Der Integralsatz von H. A. Lorentz 


Im folgenden betrachten wir ausschließlich zeitlich harmonisch-periodische 
elektromagnetische Felder, deren Zeitfaktor durch et gegeben wird. Wir nehmen 
an, daß zwei solche Felder vorhanden sind, deren komplexe Feldstärken wir bis auf 
den Zeitfaktor mit E, , H, und E,, H, bezeichnen. Setzen wir voraus, daß keine räum- 
lichen Ströme vorhanden sind, so erfüllen diese Feldstärken die Maxwellschen Glei- 


chungen 
—ix,E = rot H, ix, H — xot E, (1.1) 
wobei 
V sc Reape fan (1.2) 
c c 


Bezeichnen wir durch einen Stern den konjugiert-komplexen Wert eines Vektors, 
so folgt aus der identischen Vektorrelation 


div (E, x H; + E; x H)= H; rot E, — E, rot Hi + H, rot, E, — E;rot H, 


` mit Rücksicht auf die Maxwellschen Gleichungen (1.1) die Beziehung 


div (E, x وق‎ + E; x Hj) — 0. 


Integrieren i sie über einen durch die Fläche F begrenzten Raum R, so erhalten 
wir mit Hilfe des GauDschen Satzes den zuerst von H. A. Lorentz (1896) angege- 


benen Integralsatz 
f (OB, x HD, + (E3 x Ho f = 0 (1.3) 
EF 
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und zwar in der Gestalt, in der wir ihn weiter unten verwenden werden. Wir setzen 
dabei voraus, daß n die äußere Normale an die Begrenzungsfläche F bedeutet. Inner- 
halb des Raumes R müsen die Felder E,, H, und E,, H, im allgemeinen regulär 
sein. Gleichung (1.3) gilt aber auch, falls Sprungflächen von e und y auftreten, wenn 
nur in diesen Flächen die Tangentialkomponenten der elektromagnetischen Feld- 
stärken stetig sind. Ebenso gilt (1.3) wenn im Raume R irgend welche vollkommenen 
Leiter vorgegeben sind, da an ihren Begrenzungsflächen die Tangentialkomponenten 
der elektrischen Feldstärke verschwinden. 

Wenn das Raumgebiet R nach allen Richtungen hin sich gleichmäßig ins Un- 
endliche erstreckt, so ist es nicht notwendig in (1.3) auch den Beitrag der im Un- 
endlichen liegenden Begrenzungsfläche F zu berücksichtigen, falls nur beide elek- 
tromagnetischen Felder im Unendlichen die Sommerfeldsche Endlichkeitsbedingung 
erfüllen und ferner das eine von ihnen, etwa E,, H,, hier der Sommerfeldschen Aus- 
strahlungs- und das andere E,, H, der Einstrahlungsbedingung genügt. Überdies 
muß auch noch vorausgesetzt werden, daß die unendlich ferne Fläche F, in der 
Weise entsteht, daß eine im Endlichen liegende Fläche F gleichmäßig nach allen 
Richtungen hin ins Unendliche verlagert wird. Wir können daher F,, als eine 
Kugelfläche ansehen. 

Bezeichnet r den Radiusvektor mit einem Anfangspunkt im Endlichen, und v, 
den zugehórigen Einheitsvektor, so verlangt die Endlichkeitsbedingung daf 


lim r E — endlich, (1.4a) 
lim r H — endlich (1.4b) 


ist. Die Ausstrahlungsbedingung! fordert hingegen, daß 


lim r PEE TX n) = Oe (1.5a) 


lim r EH—r,xE)=0 5l (1.5b) 
€ 


r>00 


und die Einstrahlungsbedingung daß 


lim r PAL = Ba n) =o « (1.63) 


lim r We Yu EX 2E SS D (1.6b) 

Um zu zeigen, daß der Beitrag von مر‎ verschwindet, bemerken wir zunächst 
daß E;, H; zugleich mit E,, H, die Endlichkeitsbedingung (1.4) sowie die Einstrah- 
en Iu Lt I FETTE 


1 Bezüglich der Literaturangaben und der 


physikalischen Bedeutung der beiden Bedingungen 
(1.4) und (1.5) vgl. etwa Rubinowicz (1957). ie 
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lungsbedingung (1.6) erfüllt. Aus der Endlichkeitsbedingung (1,4a) für E, und der 
Einstrahlungsbedingung (1.6b) für E}, H; folgt sodann die Beziehung 


lim r? 2 (E, x H3) + E, x (Fo x E) | = 0 (1.7) 


und aus der Endlichkeitsbedingung (1.4a) für E; und der Ausstrahlungsbedingung 
(1.5b) für E,, H, die Beziehung 


lim ۶ v4 (Ez x Hj) — E} x (r, X E) | = 0. (1.8) 


Auf F ist nun df = r? dw, wenn dw ein Flächenelement der Einheitskugel um den 
Anfangspunkt des Vektors r bezeichnet. Da auf F, die Richtungen von T, und t 
zusammenfallen, so ergibt die skalare Multiplikation von (1.7) und (1.8) mit dem 
Vektor der Normalen m im Einklang mit unserer Behauptung das Verschwinden des 
Flächenintegrales, das in der Lorentzschen Beziehung (1.3) über die Fläche F zu 
erstrecken ist. 

In der Fassung (1.3), die wir dem Beweis des Reziprozitütstheorems zugrunde 
legen, stellt der Lorentzsche Integralsatz eine Beziehung dar, zwischen den komplexen 
Feldstirken E,, H, des einen Feldes und den konjugiert-komplexen Werten E}, H3 
der komplexen Feldstärken eines zweiten Feldes E,, H,. Man kann jedoch den Lo- 
rentzschen Integralsatz auch in der Gestalt 


JIE, x Hy), — (E, x Hy] df= 0 (1.9) 
F 


verwenden, in der direkt die beiden komplexen Feldstärken E,, H, und E,, H, auf- 
treten. Damit der Beitrag der unendlich fernen Fläche F„ zum Integralsatz (1.9) 
verschwindet, muf man jedoch fordern, daf die Felder E,, H, und E,, H, im Un- 
endlichen aufer der Endlichkeitsbedingung (1.4) beide zugleich entweder die Aus- 
strahlungs- oder die Endlichkeitsbedingung (1.5) bzw. (1,6) erfüllen. Die Fassung (1.9) 
des Lorentzschen Integralsatzes ergibt (vgl. das verallgemeinerte Reziprozitätstheorem 
für den Fall der skalaren Wellengleichung, Rubinowicz 1955) ein Reziprozitátsthe- 
orem, das ganz ebenso wie (3.16) gebaut ist, in dem jedoch in den Produkten der 
Koeffizienten A und a,, die Indizes m entgegengesetzte Vorzeichen haben. 


8 2. Entwicklung eines elektromagnetischen Feldes nach Multipolfeldern 


۱ 2 
Verwenden wir sphärische Polarkoordinaten 7, A. p und setzen k? = x, مر‎ 


so stellen nach P. Debye (1908, 1909) E", H” und E"*, H 


Kr = + kêr IIe, Ep = 0, 
i 
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20 HOSTE! m 1 ۶ 
o r 09 0r ? زو‎ gs *rsind 99 " 
el mg 
el __ i Air ft à E ees me UU 
۶ rend 7 03 
(2.1a) (2.1b) 
H0 pre LO 
r — Va r or? V 
1 or fiel 1 92 r ۶ 
ee er 
He = "e ind 9g ' ds r ddr ' 
ge 4 1 ar ۷ á: 1 ۶ 
PNA AT DD ? rsin dg dr 


zwei zeitlich harmonisch-periodische Lósungen der Maxwellschen Gleichungen (1.1) 
dar, falls die Potentiale 77۶ und 7۳۳۶ der räumlichen Schwingungsgleichung ۸477 + 
+ k? II = 0 genügen. 

Diese Lósungen kónnen auch durch die Vektorrelationen 


el 
Ed = grad T7 + kar IIS, Eng = — i xur x grad IIs, 
(2.2a) (2.2b) 
: © 3 
He! = i žer x grad I”, Hre = grad PE + k? ۶ 


ausgedrückt werden. 

Wir erhalten die elektrische bzw. die magnetische Malupo a eines 
2'-Pols!, die einer auslaufenden elektromagnetischen Welle entspricht, also der 
Sommerfeldschen Ausstrahlungsbedingung (1.5) genügt, falls wir in (2.1) oder (2.2) 
II" bzw., I” gleich 


20) y, (0, g) Q3) 


setzen. Hier bedeutet Y,, (0, el die Kugelfunktion 
Vin (D, o) = ei? PP (cos 8) (2.4) 


und 7, (el eine Funktion, die sich durch die Hankelsche Funktion H® (o) in der 
nachstehenden Weise ausdrückt: 


ldY ,ee 4 
m(e) = o! = 5 4) NAI al 173 (0) (2.5) 
MELDE Ed Eye eye 
=e ? e | + ی و‎ e A |. 


Dabei ist o = kr. 


1 l= 1, 2, 3, ... entspricht somit einer Dipol-, Quadrupol-, Oktupol- u. s. f. Strahlung. 
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Eine einlaufende, also der Sommerfeldschen Einstrahlungsbedingung (1.6) 
entsprechende elektrische bzw. magnetische Multipolstrahlung eines 2/-Pols er- 
gibt sich hingegen, wenn wir in (2.3) die Funktion n, (o) durch die gewöhnlich mit 
C; (p) bezeichnete Funktion ersetzen. Sie wird für reelle g-Werte durch den konju- 


giert-komplexen Wert von 7, (o) gegeben und ist durch die Hankelsche Funktion 
Hj? (e) in der Gestalt 


1 d eh ۶ 
up = e 24) i ; -(®) HY? 3 (o) 


EERE GD ی‎ i ED 1۳۰0 - 2710 FINE 9 
Be cons ts 


darstellbar. 
Das dem Potential (2.3) mit 7, (o) oder ¢, (p) entsprechende elektromagnetische 
Feld wird im Koordinatenursprung d. h. für o = 0 unendlich. Ein elektromagneti- 


sches Feld, das die gleichen Symmetrieeigenschaften um den Koordinatenursprung 
aufweist, jedoch hier endlich bleibt, wird erhalten, wenn wir in (2.3) n, (o) durch 


% 
viel = E» Ae )۵( 


git} 1 0 9 d e HH 
13..@+) 1-3 (2) + (2) «+ 3) 


2 2 2 


ersetzen, also aus dem Potential 
"7 Yin (9, 9). 


bezeichnet. hier eine Besselsche Funktion. 1‏ )0( چبرل 
Jedes zeitlich harmonisch-periodische elektromagnetische Feld kann durch‏ 
Superposition der angeführten speziellen. elektromagnetischen Felder Gla) und‏ - - 
P (2.1b) oder (2.2a) und (2.2b) dargestellt werden. Der Konvergenzbereich einer solchen‏ 
Darstellung wird von der Lage der Quellpunkte und der Sprung- und Begrenzungs-‏ - 
flächen in dem Raume abhängen, in dem das elektromagnetische Feld gegeben ist.‏ - 
Im folgenden betrachten wir ausschließlich solche speziellen elektromagneti-‏ 
schen Felder, die nur in einem einzigen Raumpunkte Q eine durch Übereinander-‏ 
lagerung: verschiedener elektrischer und magnetischer ele entstehende Suis‏ 
lungsquelle besitzen. Bezeichnen wir mit II" (Q, P) bzw. IT (Q, Q) die ( -‏ 
lung, des zu einem solchen elektromagnetischen Felde gehörigen II* in der Um-‏ 


` gebung eines vom Quellpunkte Q verschiedenen Raumpunktes P bzw. in der Um- 


gebung von Q selbst, so erhalten wir im Falle aus dem Quellpunkt auslaufender 
: n i 
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Wellen offensichtlich 


ITA (Q, P) = > Ya 1 (Q, P) aio Yim (8, 9). (2.82) 


I=0 دیور‎ | 


174 )0, Q) = Y S (4 (Q, Q) 12 Sg + ai (Q, sui Yin, (0, 9). (28b) 


=0 m=-I 


Die Funktion 7, (o) ist nur in (2.8b) enthalten, da sie ja nur bei der Entwicklung um 
den Quellpunkt Q auftreten kann. Die Koeffizienten in den entsprechenden Ent- 
wicklungen für //"* bezeichnen wir durch einen oberen Index mg. 

Falls es sich um eine in den Quellpunkt Q einlaufende Wellenbewegung handelt, 
ist in (2.8b) die Funktion 7, (e) (2.5) duch ¢, (o) (2.6) zu ersetzen. 

Es sei hier noch hervorgehoben, daB die Achsen der Koordinatensysteme in den 
beiden Punkten Q und P keineswegs parallel zueinander angenommen werden 
müssen. 

Zusammenfassend können wir behaupten: In einem Raumgebiet, in dem 6 
und u konstant sind, läßt sich ein beliebiges zeitlich harmonisch-periodisches und 
den Maxwellschen Gleichungen (1.1) genügendes elektromagnetisches Feld in der Ge- 
stalt 


E = Eli E”, H— H" + ۶ (2.9) 


darstellen, wo E”, H*! sich nach (2.1a) aus einem Potential ZZ” und E", H” sich nach 
(2.1b) aus einem Potential ۲۲۳۶ berechnen läßt. Die Potentiale JZZ” und analog auch 
IT" können durch Reihen von der Gestalt (2.83) und (2.8b) mit 7, (e) oder C, (o) 
dargestellt werden. 


8 3. Herleitung des verallgemeinerten elektromagnetischen Reziprozitätstheorems 


Wir bezeichnen mit R ein Raumgebiet, das sich everituell auch ins Unendliche 
erstreckt und sonst von vollkommenen Leitern begrenzt wird. Die Materialkonstanten 
e und y sollen sich in R gegebenenfalls nur in Sprungflächen ändern, sonst aber im 
allgemeinen konstant sein. In R seien zwei elektromagnetische Felder E,, H, und 
E,, H, gegeben, die je eine Strahlungsquelle in dem Raumpunkte Q, bzw. Q, be- 
sitzen. Dabei soll E,, H, eine aus dem Quellpunkt Q, auslaufende und E,, H, eine 
in den Quellpunkt Q, einlaufende Wellenbewegung darstellen. Ferner setzen wir von 
diesen beiden elektromagnetischen Feldern voraus, daß in den Sprungflächen von € 
und u die Tangentialkomponenten von E und H stetig sind und daß an den aus voll- 
kommenen Leitern bestehenden Begrenzungsflüchen die Tangentialkomponenten 
von E verschwinden. Im Unendlichen sollen beide elektromagnetischen Felder die 
Sommerfeldsche Endlichkeits- und ferner E,, H, die Ausstrahlungsbedingung (1.5) 
und E,, H, die Einstrahlungsbedingung (1.6) erfüllen, 
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Im Lorentzschen Integralsatz (1.3) verschwinden dann die Integrale über die 
Sprung- und Begrenzungsflächen und auch über die „unendlich ferne Fläche“. Die 
Integration ist daher in (1.3) lediglich über zwei die Quellpunkte Q, und Q, um- 
schließende Flächen K, bzw. K, zu erstrecken, die wir als kleine Kugeln mit den 
Mittelpunkten in Q, bzw. Q, voraussetzen wollen. Es gilt also 


J GE x B3, + (E BA df = 0. (3.1) 


KıtK, 


Da nach (2.9) 
MX H,= Et x He + Ev x Hm + 0 x He’ + Br x H# (32) 


ist und ein analoger Ausdruck für E, x H, gilt, so zerfällt in (3.1) jedes der beiden 
Integrale über K, bzw. K, in vier Integrale, die wir für die Kugelfläche K, im fol- 
genden mit (I), (II), (IIT), (IV) bezeichnen. Weiter kann dann jedes dieser vier Inte- 
grale, wenn die elektromagnetischen Felder nach $ 2 in Reihen entwickelt werden, 
durch eine Doppelreihe von Integralen dargestellt werden. Wir erhalten also für das 
Integral über K, insgesamt vier Gruppen von Doppelreihen. Zuerst behandeln wir 
das Integral 


(D = f (E? x Hg), + (EF x B2, df (3.3) 


Verwenden wir räumliche Polarkoordinaten r, 9, y, deren Ursprung im Mittelpunkte 
von K, liegt, so wird (da ja n = — rjr die äußere Normale von K, als Begrenzung 
des Raumes R bedeutet): 

(Ei x Hy"), = Et, Hj, — Eio Hag. Ga 


Um zunächst ein einzelnes Glied, der (I) entsprechenden Doppelreihe zu be- 
rechnen, müssen wir voraussetzen, daß Ej, H} und Ej, H sich mit Hilfe der 


Potentiale 
II! = f, (o) Yim (0, p) bzw. IT = fa (e) Yım (0. P) (3.5) 
darstellen lassen, in denen | ۱ 
fale) = 4... 0( LO + cf, (Or, 0) Sg (3.64) 
fale) = of On 0) ZE ` (6b) 


` bedeutet. Das Potential HS wird nämlich in der Umgebung seines Quellpunktes Qi 
also nach (2.8b) entwickelt, während 112! in der Umgebung des von seinem Quell- 
punkte Q, verschiedenen Raumpunktes Q, also nach (2.8a) dargestellt werden muß. 
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Unter diesen Voraussetzungen wird gemäB (2.1), (3.4), (3.5) und (3.6) 
(E x He"), + (ET x HD, = N (Vins Yr) M (A (e), f2@) — (9-7) 


1 ët Ee ر‎ OV im OY Em 
sin? ® Ae Ae 09 930 ° 


NP Yrs) (3.82) 


MEO rek s "9? d A SE (3.8b) 


Setzen wir nun in (I) df = r? dw, wo dw das Flächenelement der Einheitskugel e 
und r den Radius der Kugel K, bedeutet, so folgt aus (3.3) und (3.7) 


EK J N (Y,, , Vin) do. (3.9) 


l,m Jm 


Um hier zunächst die Integrale über die Einheitskugel e anzugeben, schreiben 
wir die identische Relation 


div (Y; grad Y,,) = grad Y,, grad Vi + Yr, A Yim 
in sphürischen Polarkoordinaten hin. Da 


I ) 4 1 
À Yin = m 


so erhalten wir mit Rücksicht auf (3.8a) die Beziehung 
deo 9 Y im 1.29 HOTTE 
sin A 98 (sino Vim 755 25 PTT EI dp ۱ 
= N (Ym, Yim) — L (I + 1) Ym Yr. 


Bei der Integration über die Einheitskugel e verschwindet die linke Seite dieser Rela- 
tion (dw = sin à dŷ dp) und es wird daher 


LN Om Yr) do = 102-1) f Yim Yr, do = 1 (1 + 1) Än Se Kim 


wo Ky, = 2x | (Pf (cos #))? sin 9 dô (3.10) 
"$0 


das Quadrat des Normierungsnenners der Kugelflächenfunktion Yi, bezeichnet. 
Um (I) anzugeben, muß noch M (f, f;) berechnet werden. Aus (3.6a), (3.6b) 
und (3.8b) folgt 


M (f. fà = Af (Q1, Q) af (Qe, Q) M bg Jr (3.11). 


 — =» ———— 


. ۰ "m, 
setzen wir voraus, daß E, H und E>”, 
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wenn wir mit Rücksicht auf (3.10) /' = / und m’ = m setzen. Aus der Differential- 


relation für die Zylinderfunktionen 


d Ly» n 
ae =. Zn (0) + Zn-1 (0) 


sowie der Funktionalgleichung 


Ltd H9 (9) — J, (9 H0, (9) = = 


folgt nun die Beziehung 
d dai. i. 
N do yes i do US ceto 


Es wird somit 


m (2, 2) - - tr 
0 0 0 


(3.12) 


Wir erhalten daher aus (1.2), (3.9), (3.10), (3.11) und (3.12) für (T) (3.3) den Ausdruck 


oo l 
(D = — > y ]0 + yz 44 (Q1, Q1) af (Q2, Q1) Kim: 


` 1=0 m=- 
Auf dem gleichen Wege ergibt sich für das Integral (vgl. (3.3)) 
un = f (BF x Hy), + (B3 x DU df 
Kı 


die Reihenentwicklung 


oo 1 
m=- Y Y 101 1(1/ © 4% (Q1: QD af (a O) Kan: 


I=0 m=-I 


Hingegen kann mafı zeigen, daß die beiden Integrale 


am = f (Er x By), + (BFF x Hi) d und 
Ki 


qv) = f (E x Hi), + (EY x BM] df 
Kı 


(3.13) 


(3.14) 


verschwinden. Wir führen den Nachweis für das Integral (III) durch. Es genügt ein 
einzelnes Glied der Doppelreihe zu betrachten, die (III) darstellt. Zu diesem Zwecke 


In = fs (e) Yim (0: g) bzw. II? = fa (0) Yr (9, 9) 


Hy’ durch die beiden Potentiale 
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gegeben werden. Man erhält dann 


(ES x Hz*), + (BFF x Hj), = O (Yu. Yrm) P (fafa) 


1 ( Yeo Yi, ON SA 


O Qin, Tu) = re Ee 


Pang en‏ از 
P (USE E DE "ET‏ 


Um zu zeigen, daß (IIT) verschwindet, genügt es den Nachweis zu erbringen, daß 


10 s Yew) O 


ist für beliebige l, m und J’, m’. Dies folgt aus der Tatsache, daß wegen rot grad = 0 
die identische Beziehung 


rot (Y,, grad Fo) = grad Y,, x grad Ys? 


gilt, deren r-Komponente in sphärischen Polarkoordinaten ausgedrückt 
ad 1 * 
rot, Bän grad Yim) SC 7 2 4( ) Yrm) 


lautet. Damit ist unser Ziel erreicht. Es ist ja dann 


| Ota, Yin) do = [ rot, (Ym grad Yî) d 
& K, 
und das hier auftretende Integral über K, verschwindet wegen div rot = 0. 
Ebenso kann man zeigen, daß auch (IV) = 0 ist. 
Im ganzen erhalten wir somit mit Rücksicht auf (3.13) und (3.14) für den Beitrag 
der Kugel K, zum Integral (3.1) den Ausdruck 


o 1 
mw 3 ES 1 0 + 1) Ks, T (Q1, Q1) at, (Qa, O) + 
1 


I=0 m=— 


+ Vi ARE (O1. QD af" (Qs, oj]. (8.15) 


Der Beitrag der Kugel K, zum Integral (3.1) berechnet sich ganz ebenso wie 
der Beitrag der Kugel K,. Dabei ist nur zu berücksichtigen daß das elektromagne- 
tische Feld Ej, H, durch eine im Quellpunkt Qs einlaufende Welle gegeben wird 
und daher in seinem Potential die Funktion ¢, (e) (2.6) an Stelle von 7, (o) (2.5) in 
E,, H, tritt. Die Rechnung ergibt für den Beitrag der Kugel K, zum Integral (3.1) - 
den konjugiert-komplexen Wert des Ausdruckes (3.15) in dem die Indizes 1 und 2 
vertauscht und das Vorzeichen geändert wurde. Man erhält daher aus (3.1) das ge- 
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în 


suchte elektromagnetische Reziprozitätstheorem in der Gestalt 


eo 1 x 
Y Y ung, Wi Ad (Qu, Qs) af (Qas Q) + 


1-0 m=-—I 


a y ABE (Qu 02 که‎ (0, 00] = (019 
۲ 1 


Li‏ مه 
cox Der UE Aih (Qa. Qa) aff, )0 ON +‏ 
:2 


I=0 m=— 
+ 9/8 am (On O9 aff Ov o»). 
2 


Der Tatsache, daß & und y in den beiden Raumpunkten Q, und. Q, verschieden sein 
können, wurde in (3.16) dadurch Rechnung getragen, daß diese beiden Material- 
konstanten mit den entsprechenden Indizes versehen wurden. ` . 


8 4. Diskussion der Rekursionsformel 


Für die folgenden Überlegungen setzen wir voraus, daß das elektromagnetische 

Feld E,, H,, das eine in ihrem Quellpunkte Q, einlaufende elektromagnetische _ 
"Wellenbewegung darstellt, nur eine einzige, einem speziellen Indexpaar Iq, Mo ent- 
sprechende elektrische oder magnetische Multipolquelle mit der Amplitude Att (Qs O5) 
bzw. ATE, (Qa, Qe) besitzt. Im ersten Falle verschwinden dann alle übrigen andere 
. Indizes 1, m tragende Af, (Qe, Qa) sowie ausnahmslos alle 477 (Q,, Qa). Im zweiten 
x Falle sind alle mit anderen Indizes versehenen AS (Q2, Q3) = 0, sowie ausnahmslos 


3 A ETATE gode 


P 


4 
NGA win PUR A8 
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Es muß darauf hingewiesen werden, daß die Koeffizienten afl, (Qa, Q1) und 
a? (Q,, Q,) i jeder der beiden Relationen (4.1) und (4.2) zu elektromagnetischen 
Feldern mit verschiedenen Strahlungsquellen gehóren und daher auch verschiedene 
Werte haben. Hingegen kann sowohl (4.1) als auch (4.2) selbstverstándlich auf das 
gleiche elektromagnetische Feld Ej, H, angewendet werden. 

Um zu zeigen, in welcher Weise das Reziprozitütstheorem (4.1) bzw. (4.2) zur 
Berechnung von Multipolfeldern verwendet werden kann, setzen wir ferner voraus, 
daß das elektromagnetische Feld E,, H,, das in (4.1) bzw. in (4.2) auftritt, bekannt ist. 
Wir können dann auch die Koeffizienten af, (Qa, Q,) und of? (Qs, Qı), die die Ent- 
wicklung dieser Felder in dem Raumpunkte Q, darstellen, sowohl im Falle (4.1) als 
auch im Falle (4.2) als bekannt ansehen. Denken wir uns nun ein elektromagnetisches 
Feld E,, H, durch die Amplituden A3. (Q1, Q1), AR (Oh Q1) seiner Multipolquellen 
im Quellpunkte Q, definiert. Wir können dann mit Hilfe von (4.1) bzw. (4.2) die 
Koeffizienten at ns (Qi. Qe) bzw. az£ (Q1, Q2) seiner Entwicklung nach elektrischen 
bzw. magnetischen Multipolfeldern im Raumpunkte Q, angeben. 

Will man nun den Verlauf des Feldes E,, H, in der Umgebung des Punktes Q, 
bestimmen, so muß man alle Koeffizienten a^ (Q,, Q) und agë (Q1, Q,) der Ent- 
wicklung von E,, H, im Punkte Q, kennen. Es ist dann erforderlich die beiden Rezi- 
prozitätstheoreme (4.1) und (4.2) für alle Jj, mg anzuwenden. Man muß dann also 
die Felder E,, H, mit elektrischen bzw. magnetischen Multipolen beliebig hoher 
Ordnung zur Verfügung haben. 

Viel einfacher ist das Problem die elektrischen und magnetischen Feldstärken 
des Feldes E,, H,, das beliebig vorgegebene Multipolquellen im Punkte Q, besitzt, 
nur in dem Punkte Q, zu berechnen. Diese Feldstärken werden, wie in 85 klar werden 
wird (vgl. (5.8)), durch die Koeffizienten af, (Q1, Q2) und af (Q1, Qa) (m = 0, + 1) 
bestimmt. Wie aus (4.1) und (4.2) hervorgeht, ist es zu diesem Zwecke nur erforder- 
licn Felder E,, H, mit nichtverschwindenden 47, (Q2, Qa) und A (Qz, Qs) 
(m = 0, + 1) zur Verfügung zu haben d. h. elektromagnetische Felder, die einer 
elektrischen und einer magnetischen Dipolstrahlung entsprechen, ohne Rücksicht 
auf die Ordnung der elektrischen und magnetischen Multipolquellen des Feldes 
E,, H,, das berechnet werden soll. 


8 5. Berechnung einiger Koeffizienten af, und of? aus den elektromagnetischen 


F: eldstärken 


Bei den Kosselschen Interferenzversuchen beobachtet man praktisch im Un- 
endlichen die Beugungsfransen, die durch die Beugung der von einem Atom eines 
Kristallgitters emittierten Röntgenstrahlung an anderen Atomen des gleichen Kristall- 
gitters entstehen. Im Falle, wo die Röntgenstrahlung von einem elektrischen Dipol 
herrührt, hat Max von Laue (1935, 1937, 1941, 1948) eine sehr elegante Lösung dieses 
Problems gegeben. Er hat nämlich auf die bekannte Lösung des Problems der Re- 
flexion einer ebenen Röntgenwelle an einem Kristallgitter das Reziprozitätstheorem 
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für die elektrische Dipolstrahlung angewendet. Da man annehmen kann, daß die 
einfallende ebene Welle durch eine im Unendlichen liegende Dipolquelle verursacht 
wird, kann man das elektromagnetische Feld des Reflexionsproblems als das Feld 
ansehen, das in unserer Formulierung des verallgemeinerten Reziprozitätstheorems 
mit E,, H, bezeichnet wird. Durch Anwendung des Reziprozitätstheorems für die 
elektrische Dipolstrahlung kann man sodann das elektromagnetische Feld E,, H, 
einer in einem Gitterpunkte liegenden Dipolquelle am Orte des Quellpunktes von 
E,, H, d. h. im Unendlichen angeben. 

In der gleichen Weise kann man unter Zuhilfenahme der gleichen Lösung des 
Reflexionsproblems durch Anwendung des für beliebige Multipolquellen verallge- 
meinerten Reziprozitätstheorems den Fall erledigen, wo ein dem Kristallgitter ange- 
hörendes Atom irgend eine Multipolstrahlung aussendet. In einem jeden solchen 
Falle steht man vor der Aufgabe die im betrachteten Multipolfalle im Reziprozitäts- 
theorem (4.1) oder (4.2) auftretenden Koeffizienten aj. (Qs, Q1) und aff (Qs, Q1) 
aus den Feldstärken E,, H, zu bestimmen. 

Im Falle der elektrischen und magnetischen Dipol- und der elektrischen Qua- 
drupolstrahlung, die bei vielen physikalischen und technischen Problemen auftreten, 
ist die Kenntnis der Koeffizienten 


afl, (Qa, Q) ) — 1, 2, Gl) a (Qu ON  (5-1b) 


erforderlich. 

Die Koeffizienten af (Q,, Q1), aod (Qe, Q1) können wir gleich Null setzen, da 
sie ebenso wie die Koeffizienten 47^ (Q;, 0), A (Qj, Q) @ = 1, 2) in der Entwick- 
lung von E; H; (i = 1,2) überhaupt nicht auftreten. In den Entwicklungen der zu- 
gehörigen Potentiale enthalten die entsprechenden Glieder mit Rücksicht auf Yoo (Ô, p) 
— const nur die Funktionen y, (o), no (e) oder Co (o). Gemäß (2.1) verschwinden 
daher die Radialkomponenten E? und Hy", weil ja Zl=0 oder H7*—0 
geradezu die Differentialgleichung darstellt, der diese Funktionen genügen. Die übri- 
gen Komponenten von E", H” und E", H"* sind aber gemäß (2.1) gleich Null, 
weil die Glieder mit den Koeffizienten a, und ag mit Rücksicht auf Yoo (8,9) — 
const von 9 und g nicht abhängen. Diese Glieder besitzen daher eine Bedeutung nur 
für die Entwicklung der Potentiale IT und IT”, nicht aber für die Entwicklung des 
elektromagnetischen Feldes. ; 

Der Zahlenwert der Koeffizienten (5.1) hängt selbstverstándlich von der Nor- 
mierung der in Y,„ (0, p) (2.4) enthaltenen zugeordneten Kugelfunktionen Pr (EI ab. 
Wir entscheiden uns für die von C. G. Darwin angegebene Normierung, setzen also 


pii; Pi p Pi 9%; PO 3 بقع‎ Pj = (1— E9538 P1730 — E) 
i ۱ > eg put (5.2) 
und nehmen dementsprechend an, dab p^" (8) = (SDE Em (£) ist, 
Um zunächst die „elektrischen“ Koeffizienten (5.1a) anzugeben, ist es notwendig 
die Feldstärken ES, H? nach den Koordinaten x, y, 2 bis einschließlich der linearen 
Glieder zu entwickeln. Aus (2.2a) ist ersichtlich, daß hierzu nur die Kenntnis der 
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linearen und der quadratischen Glieder von JI” erforderlich ist. Man muß daher in 
der Entwicklung (2.82) des Potentials II* in der Umgebung eines regulären Punktes 
in den Funktionen y, (p) nur noch dritte Potenzen von o berüchsichtigen. Gemäß 
(2.7) genügt es daher y, (o) = 07/3, Ya (o) = o?/15 und alle übrigen y, (o) = 0 zu 


setzen. Führen wir die Bezeichnung 


ON ey Vin (9, 9) (5.3) 


ein, so haben mit Rücksicht auf (5.2) die von uns benötigten 7774, wenn wir sie bis 
einschlieBlich der quadratischen Glieder in x, y, z angeben, die Gestalt 


IE, = — ad Qik(x— Ly), IE toos IT], = af, 3 k ) + iy), 


3x d$ a 3 käte-— iy), Lie d$ 5 k? (x — iy) z,‏ = و11 
(x + ois‏ چ af‏ = رو ,)9 — 2 — 222( 2 وه = IL‏ 
als ty)? | (5.4)‏ و11 


Andere IZ kommen nicht in Frage, da ja das م7‎ entsprechende elektromagne- 
tische Feld verschwindet und die HS, (5.3) mit IS 3 vom mindestens dritten Grade 
Eh 2 Sind: 

Um die Formeln in (5.4) nicht allzu lang werden zu lassen, haben wir bei den 
hier auftretenden Koeffizienten af, (Us, Qı) die Klammerbezeichnungen (Q,, Q;) 
fortgelassen. Dies soll aus dem gleichen Grunde auch bei allen weiter unten anzu- 
gebenden Formeln geschehen. Man vergesse also nicht, daß sie sich auf ein elektro- 
magnetisches Feld mit dem Quellpunkt in Q, beziehen, das wir in dem Raumpunkte 
Q, oder in seiner Ungebung betrachten. 

Bezeichnet man die IZ, (5.3) gemäß (2.2a) entsprechenden Feldstärken mit 
Hi, so erhält man aus den Potentialen (5.4) für sie die Ausdrücke 


E? 


in 
Ej Recte Gras 5» (E— ij). Hj, u des ik ee ij 

| + k (x — iy)], 
Ej, = وه‎ ikk, Hy, = af) 5k x, [iiy — jix, 
Ei = ad $k (i + i), Hi = at 3 k x, [(È + ij) z — 

— k (x + iy). 

Ee $k (i — ij) (x — iy)  Hł_,=0, ۱ (5.5) 
E ee Sé a2 و‎ E? ]6 — i$) z + a, 


o ex zd nts 
Es, = an: k? [—ix— jy + k22], He, = 0, 
Ej =o%,2[@+ij)z+k(x+iy)], Hf =0, 
Ez =a tk 06+ ( (x +i) Hf =0, 


falls wir uns auf die in x, y, z linearen Glieder beschränken. 
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Um auch die magnetischen Koeffizienten a%€ (Q,, Q,) (5.1b) zu erhalten, müßten 
wir das magnetische Potential ۲۳۳۶ in Potentiale [77# aufteilen, die den elektrischen 
Potentialen IZ, (5.3) entsprechen und sodann die aus den IT}: gemäß (2.2b) abzu- 
leitenden Feldstärken EF, H7 berechnen. Es ist jedoch gar nicht notwendig diese 
Rechnungen durchzuführen, da, wie der Vergleich von (2.2b) mit (2.2a) lehrt, wir 
diese Feldstärken erhalten, wenn wir an Stelle von Ei}, bzw. Hr, die Feldstärken H7# 
bzw. — Ej treten lassen. Dabei ist gleichzeitig x, durch x, zu ersetzen und man hat 
statt af! die Koeffizienten a7 (Q,, O zu verwenden. Für den Teil II der vorliegenden 
Arbeit benötigen wir von den elektrischen Feldstärken E7* nur solche mit den in 
x, y, z linearen Gliedern und von den magnetischen Feldstärken Hf nur solche mit 
konstanten Gliedern. Die in dieser Näherung nicht verschwindenden Feldstärken 
Enz, Hy? ergeben sich aus (5.5) zu 


me = — at ikx,[ü—-ij)z—k&-iy], MT, atb — OH, 
Ej = — aff 5k %, (iiy — jix), Hi: - ۶ (5.6) 
Em ——oWlkx[ü-ijz-k(c-iy)h Hi -etgk-cij. 


Das gesamte vom Quellpunkte Q, herrührende elektromagnetische Feld erhalten 
wir in dem Raumpunkte Q,, wenn wir die Feldstärken (5.5) und (5.6) zueinander 
addieren: 

E k? (x — iy) — ag 4 ` (ër‏ و Zr sh + af‏ تچ سا سل 
gy à À? (x + iy) — am 3 2 — avg Ek x, (x + iy),‏ — 
k? (x + iy) +‏ £ وه + 2 ig at Eg at ik = ab S k? (x — iy) + ag, $k‏ 
amet k x, (x — iy) — att S ko, z,‏ + 
k” (x + iy) +‏ 5 ره + E,— ad 2h — ad 48 k? (x — iy) + ago 5 k?z‏ 
ars lku (x — iy) + att 5 x, (x + (۰ |‏ + 
mi ET e 2‏ ` . 
H,— iH, = Oe km; = aM ik (57‏ ماج 0۳۴ - H,--iH,—‏ 
Wir haben dabei bei der elektrischen Feldstärke E nur die in x, y, 2 linearen und bei‏ 
der magnetischen Feldstärke H nur die konstanten Glieder berücksichtigt.‏ 
Aus (5.7) erhält man entweder unmittelbar oder durch Nullsetzen von x, y, 2‏ 


0 N e Mgr eae 
oder auch durch Bildung der Ableitungen = +i dy? de fiir die Koeffizienten (5.1) 


die Ausdriicke 
p 8otk= ES 404 k= bn if, 
$a] k = E, + iE, 2 موه‎ = E, و‎ fk = E, y? 


A 4 ng h El 
- ۵ وا دی‎ H, + iH,» ga b Hs şaf k= H, — iH, 
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9 . 9 e 12 d d 
bei a$ ak = (oz +i HI (E dE. ou per k3 = CS (E,+iE,) + 


d ) 4 IE, 12 d 
+ (£i) = af, kt = —* g ay B= (E iE) + 


d 123: 24. a 0 .9 
SR ۳ =] iy) Be FX 2 — EI : با‎ iE. (5.8) 


Die Koeffizienten (5.1) stellen das durch die Quelle Q, in dem Raumpunkte Qj 
verursachte Feld dar. Die in (5.8) auftretenden elektromagnetischen Feldstärken 
und ihre Ableitungen entsprechen daher dem Felde E,, H, in dem Raumpunkte Q,. 
Selbstverständlich kann man die Beziehungen (5.8) auch dazu verwenden um die 
Koeffizienten aj. (Oh, Q) 0 = 1, 2) und af (Q,, Qa) des Feldes E,, H, mit dem 
Quellpunkt in Q, in dem Raumpunkte Q, darzustellen. 

Da zwischen den elektromagnetischen Feldstärken und ihren Ableitungen die 
durch die Maxwellschen Gleichungen (1.1) geforderten Beziehungen bestehen, kann 
man die Koeffizienten (5.1) auch noch in anderer Weise durch die elektromagneti- 
schen Feldstärken ausdrücken, als dies in (5.8) geschehen ist. Wir haben uns für die 
einfachste Darstellung entschieden, in der die „elektrischen“ Koeffizienten ati durch 
die elektrischen und die „magnetischen“ Koeffiziehten aps durch die magnetischen 
Feldstärken gegeben werden. 

Bemerkenswert ist vielleicht eine Ähnlichkeit im Aufbau der Koeffizienten af}, 
und az (5.8) und im Aufbau der konjugiert-komplexen Werte der räumlichen Kugel- 
funktionen, die die gleichen Indizes Z, m tragen und die z. B. in den Potentialen 777 
(5.4) auftreten. Die Koeffizienten af, bzw. of? werden nämlich aus den entsprechen- 
den konjugiert-komplexen ráumlichen Kugelfunktionen erhalten, wenn man in ihnen 
die Koordinaten x, y, z teils durch die Feldstürken E,, E,, E, bzw. Ho H,, H. teils 
durch die Differentialoperatoren 9/9x, 9/9y, 9/9z ersetzt. Dies gilt auch im Falle des 
Koeffizienten v. wenn man beachtet, daß man mit Rücksicht auf div E = 0 ihn 
auch in der Gestalt 


darstellen kann. Um jedoch festzustellen, ob dieser Analogie eine tiefere Bedeutung 
zugrunde liegt, müßte man die den höheren Multipolen entsprechenden Koeffizienten 


ag und az durch die elektromagnetischen Feldstärken und ihre Ableitungen aus- 
drücken. 
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SUR LES SPECTRES RAMANIENS DES SOLUTIONS DE PYRI- 
DINE ET DE 2,06 — LUTIDINE DANS LES ACIDES SATURÉS. 


Par ROMAN MIERZECKI 


Institut de Physique Expérimentale de l'Université de Varsovie 


Laboratoire d'Optique Moléculaire de l'Institut de Physique de l'Académie Polonaise des Sciences 
(Reçu le 19 Mai, 1959) 


On a étudié les spectres ramaniens des solutions de différentes concentrations de la 
pyridine dans l'acide acétique et l'acide formique, ainsi que les solutions équimolaires 
de la pyridine dans l'acide butyrique et de la 2,6 — lutidine dans les mémes acides saturés. 
On a constaté, que les variations du spectre de la pyridine et de la 2,6 — lutidine sous l'in- 
fluence des acides formique et butyrique sont analogiques aux variations observées 
auparavant dans les spectres des solutions de la pyridine dans l'acide acétique. (Mie- 
rzecki 1953, 1955a). En particulier on observe des raies nouvelles correspondant aux 
pulsations du noyau (1007 K environ) et aux oscillations des liaisons C—H (3070 K en- 
viron) dans les molécules de la pyridine qui forment des complexes avec les molécules 
des acides. 

Les mesures d'intensité de la raie 1003 K en fonction de la concentration des solu- 
tions indiquent, que le complexe qui se forme dans les solutions étudiées se compose 
d'une molécule de pyridine et d'une molécule d'acide. 

On a montré aussi, que l'apparition de la raie de fréquence 1020 K qui a été observée 
dans les spectres de la pyridine dans l'acide acétique (Mierzecki 1955a, 1956) dépend 
de la présence d'une certaine quantité d'eau dans les solutions étudiées. Cette raie forme 


une modification d'une des raies de la pyridine. 


I. Introduction 


Les systémes bases pyridiniques — acides saturés présentent le cas, où il ge 
des interactions intermoléculaires prédominantes. Dans des systèmes pareil méme 
la formation de complexes est possible. En vue de ces effets plusieures chercheurs 
ont étudié dans le systeme pyridine — acide acétique les changements de différentes 
propriétés physico-chimiques en fonction de la concentration. A parür des points 


e ک۹‎ 
۲ A ; c nflec- 
particuliers des courbes qui représentent ces fonctions (extremums, points d'i 


tion), ils ont essayé de déduire des conclusions concernant la composition des com- 


plexes. Selon la propriété étudiée, différents auteurs proposent différentes composi- 


(41) 
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tions des complexes. La proposition la plus souvent rencontrée c'est 1 molécule de 
pyridine pour 4 molécules d'acide. D'aprés d'autres propositions les nombres de 
molécules de ces deux. substances dans un complexe se comportent comme 3: 2, 
2:3, 1:3. Les études crioscopiques seules (Swearingen, Randal, Ross 1934) mon- 
trent, que le complexe qui se forme dans les solutions de la pyridine et de l'acide 
acétique, doit étre composé d'un nombre égal de molécules de pyridine et d'acide 
acétique. 

Les fortes interactions intermoléculaires qui existent dans les solutions doivent 
se manifester aussi par des changements dans les spectres ramaniens. Les résultats 
des mesures de fréquence des raies dans les spectres des solutions de la pyridine et 
de l'acide acétique de différentes concentrations obtenus par l'auteur ont été publiés 
auparavant (Mierzecki 1953, 1955 a, b). Les résultats obtenus en méme temps par 
Lakshmanan (1954), Chiorboli (1957) et Schwab et Glatzer (1957) sont en accord 
avec les résultats de l'auteur. 

Il était intéressant d'étudier aussi les spectres des solutions de la pyridine et 
de la 2,6-dimethylpyridine (2,6-lutidine) dans d'autres acides saturés pour chercher 
des changements analogues des spectres. On a étudié alors les spectres des solutions 
équimolaires de la pyridine dans les acides formique et n-butyrique, ainsi que les 
spectres des solutions équimolaires de la 2,6-lutidine dans les acides formique, acé- 
tique et n-butyrique. 


II. Expérimentale 


a) Mesures de fréquence des raies 

Les spectres ont été étudiés par la méthode photographique. Les détails de pré- 
paration des solutions de la pyridine et de l'acide acétique, ainsi que de mesure de la 
fréquence des raies ont été décrites auparavant (Mierzecki 1953, 1955a, 1957). Les 
solutions équimolaires des bases pyridiniques avec des acides ont été préparées en 
ajoutant á la pyridine distillée une quantité équivalente d'acide. 

On a étudié aussi les spectres des solutions de la pyridine et de l'acide formique 
contenant 7, 9, 20 et 35% mole de pyridine. Les substances ont été soigneusement 
deshydratées et des précautions spéciales ont été prises pour éviter l'absorption de 
l'eau de l'atmosphére. 

b) Mesures d'intensité des raies 

On a étudié la variation d'intensité de quelques raies dans les E des solu- 
tions de la pyridine dans l'acide acétique en fonction de la concentration. 

Les méthodes photographiques de mesure et de comparaison de l'intensité des 
raies Raman dans les spectres des solutions de différentes concentrations sont trés peu 
précises. Dans les études présentées il a été impossible d'utiliser la méthode d'étalon 
intérieur, car l'addition d'une troisiéme substance pourrait influencer les interactions 
interinoléculaires. De méme la comparaison avec l'intensité d'une raie de mercure 
n'est pas satisfaisante, car dans les spectres des solutions de différentes concentra- 


LÀ 
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tions son intensité n'est pas constante. Cette observation faite au cours de ce travail 
n'est pas étonnante, car l'intensité des raies de Rayleigh dépend de plusieurs facteurs 
qui n'ont pas la méme valeur dans toutes les solutions étudiées, Pour cette raison on 
a pu suivre seulement les changements d'intensité des raies les plus intenses en fonc- 
tion de la concentration. En particulier on a étudié l'intensité des raies de fréquences 
989 et 1029 K de la pyridine pure et celle de la raie nouvelle de fréquence 1003 K. 


Ces trois raies forment une bande commune. 


Les spectres ont été photographiés à l'aide d'un spectrographe Huet B II. La 
concentration des solutions a été la méme, qu'auparavant à savoir: 100, 91,5, 80,8, 
73,2, 61,0, 48,3, 43,3, 39,8, 37,8, 30,4, 22,3, 21,6, 14,6 et 6,5%, mole de pyridine. 
La fente du spectrographe a été élargie jusqu'à 3,5 K. À cette largeur de la fente on 
a obtenu des raies suffisamment fortes et en méme temps suffisamment fines pour 
séparer les trois raies formant la bande. À l'aide d'un microphotomètre enregistreur 


de Moll on a obtenu des courbes photométriques. Pour séparer les trois raies on a re- 


calculé ces courbes en diagrammes d'intensité à l'aide des courbes caractéristiques 
pour chaque plaque. En supposant, que les raies composantes ont une forme 
symétrique on a pu facilement les tracer. Elles sont indiquées sur la fig. 1 par 
pointillés. Les aires limitées par ces courbes sont proportionelles à l'intensité 
intégrale. 
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t dans la région de 1000 K dans les spectres des solutions de la 
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pyridine et de l'acide acétique de différentes concentrations d' 
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III. Résultats 


Les résultats des mesures de fréquence des raies dans les spectres des solutions 
équimolaires de la pyridine sont indiqués au tableau 1. Ces résultats sont représentés 
également sous la forme d’un diagramme (fig. 2). 


Tableau I 


Fréquences des raies dans les spectres des solutions équimolaires de la pyridine (en K) 


solvent 


pyridine pure* 


601 w 


650 s 


888 w 
989 vs 
1029 vs 
1215 m 
1573 m 
1580 m 
1596 w 
3055 vs 


solvent © 


7000 


CH,COOH* 


7500 


n-C,H,COOH 


602 wv 


650 m 
861 w 
887 w 
989 w 
1003 s 
1029 m 
1214 w 
1573 w 
1596 w 


3067 m 


2000 3000 K 


Fig. 2. Diagrammes des spectres des solutions équimolaires de pyridine. 


* D'aprés Mierzecki (1955 a, b) 
(vs — trés forte, s — forte, m — moyenne, w — faible, vw — trés faible). 
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Le tableau 2 compare les fréquences des raies Raman de la 2,6-lutidine obtenue 
par différents auteurs. Le tableau 3 et la figure 3 montrent les fréquences des raies 
dans les spectres des solutions équimolaires de la 2,6-lutidine. 


Les résultats des mesures de fréquence des raies dans les solutions deshydratées 
de la pyridine dans l'acide formique sont indiqués au tableau 4. 

Les valeurs relatives de l'intensité intégrale des raies 989, 1003 et 1029 K dans 
les spectres des solutions de la pyridine et de l'acide acétique d'aprés trois plaques 
sont représentées sur la figure 4 a b c. Pour toutes les plaques on a admis l'intensité 
de la raie nouvelle dans la solution contenant 48%, mole de pyridine égale à 100. 
Les résultats sont précis à 25% prés. 

On a introduit des corrections pour l'influence du changement de l'indice de 
réfraction dans les solutions de différentes concentration sur l'intensité mesurée. 
Suivant Woodward et George (1951) on a multiplié l'intensité mesurée par l'indice 
de réfraction au carré. En se basant sur les valeurs de l'indice de réfraction pour les 


Tableau II 


Spectres de la 2,6-lutidine d'aprés différents auteurs 


Manzoni Ansidei Herz, Kahovec E Travail 
(1940) Kohlrausch (1943) présent 


197 m 
287 m 
427 vw 
536 m 
553 m 
714 s 

892 vw 
994 vs 
1038 vw 
1092 w 

1155 vw 
1262 m 
1274 w 

1373 m 
1404 vw 
1442 vw 
1578 m 
1604 vw 
2922 s 


3061 m 
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solutions de la pyridine dans l’acide acétique d’après Pushin et Matavuyl (1932) on 
a construit un diagramme qui représente la variation de l’indice de réfraction en 
fonction de la concentration. De ce diagramme et du diagramme représentant les 


Tableau III 
Fréquences des raies dans les spectres des solutions équimolaires de 2,6-lutidine 


solvent 


CH4COOH n-C,H,COOH 


lutidine pure 


197 m 201 m 201 m 
287 w 288 w 291 w 
427 vw — == 
536 m 37 DT 538 m 
553 m 554 m = 
— 561 w 561 m 
714 s 717 m 719 m 
= 874 w 864 w 
892 vw — — 
994 vs 994 w 995 vw 
— 1008 s 1008 m 
1038 vw LI — 
1092 w 1097 vw 1098 vw 
1155 vw 1157 vw 1160 vw 
1262 m 1266 m 1268 m 
1274 w dm — 
1373 m 1378 m 1380 w 
1404 vw = — 
1442 vw == 1441 vw 1450 vw 
1578 m 1578 w 1579 w 1580 m 
1604 vw 1597 vw 1597 vw 1593 m 
2922 s 2927 s 2926 s 2935 s 
3061 m 3066 w 3065 m 3068 w 
solvent 0 500 1000 1500 2000 3000 K 


Fig. 3. Diagrammes des spectres des solutions équimolaires de 2,6-lutidine. 
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Tableau IV 


Spectres ramaniens des solutions de la pyridine et de l'acide formique à concentration de pyridine 
inférieure à l'équimolaire 


acide 


7% mole 9%, mole 


formique* 20% mole 35% mole pure** 
ENETLIAES Em Fk Ps 
221 vw 221 vw 5 
= 611 w 611 w 610 w = 601 w 
== 639 w 639 w 639 w 642 s — 
= = = — 651 m 650 s 
679 m 680 m 680 m 680 m lt och 
695 w 701 m 703 m 701 m 696 w as 
e = = — 994 vw 989 vs 
= 1011 vs 1011 vs 1011 vs 1010 vs — 
== =; = — 1022 m = 
= 1029 w 1029 w 1028 w 1033 sh 1029 vs 
1048 w 1061 w 1061 w 1062 w 1063 w == 
1205 m 1205 m 1205 m 1205 m 1208 m — 
= == = 1253 w 1254 w 1215 m 
1345 vw 1345 w 1348 w 1349 w 1347 w — 
1395 s 1398 s 1398 s 1399 s = = 
— — — — 1578 w 1573 m 
SÉ = 2 == == 1580 m 
— —. 1604 vw 1604 vw 1604 w 1596 w 
— — 1640 m 1640 m 1636 m — 
1670 s 1672 s 1680 s 1680s — = 
Br EZ 1716 s 1719 s 1714 w — 
1725 m 1734 s 1730 s — 1740 w = 
2961 vs 2956 vs 2952 vs 2939 vs 2912 s = 
E =, = = 3072 vs 3055 s 
— 3105 vw 3100 vw 3091 vw 3094 m = 
- + = ۱ 3107 ۷ — — 


valeurs de l'intensité intégrale mesurée (fig. 4 b) on a pris les valeurs pour les solu- 
tions contenant 100, 90, 80, 70, 65, 60, 55, 50, 45, 40, 30, 20, et 10% mole de pyri- 
dine et on a calculé les valeurs corrigées de l'intensité intégrale (fig. 5). 


I V. . Discussion 


Dans les spectres ramaniens des solutions de la pyridine et de l'acide acétique 


on observe des changements de fréquence qui dépendent de la concentration de la 
pyridine, ainsi que l'apparition de raies nouvelles absentes dans les spectres des sub- 
stances pures. On observe en particulier des raies nouvelles de fréquence 628, 1003, 


et 3068,5 K parallélement à la disparition partielle des raies 604, 989 et 3055 K de la 


* D'aprés Chiorboli (1957). 
** D'aprés Mierzecki (19552). 
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pyridine pure. Ces derniéres correspondent alternativement aux oscillations 6 a, 
1 et 2 (Kline, Turkevich 1944, Corrins et coll. 1953) (fig. 6). On observe aussi des 
changements dans la région 1600 K. Dans les spectres des solutions équimolaires 
de la pyridine dans d'autres acides saturés on observe des changements analogues 
(Lakshmanan 1954, Mierzecki 1957, Chiorboli 1957, Schwab et Glatzer 1957). 

a) Changements dans la région de 1000 K. 

L'apparition d'une raie nouvelle dans cette région a été observée dans les spectres 
des solutions de la pyridine dans différents solvents. Les fréquences de ces raies pour 


les solutions équimolaires sont indiquées au tableau 5. — ` 


7 
U 
e plaque N°55 
200 o plaque N° 58 
x plaque N* 61 
X 
150| i j zx 


Ww 8m 
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150 e plague NS 55 
o plaque N° 58 
x plague N° 61 


100 


I00 80 60 4D 20 0 
% mol ۷ 


Fig. 4b 


Différents auteurs ont observé dans les spectres des solutions aqueuses de la 
pyridine la raie de fréquence 998 K au lieu de la raie 989 K de la pyridine pure. Les 
changements des courbes photométriques pour les spectres des solutions aqueuses 
contenant 79 et 66,59/, mole de pyridine (fig, 7) montrent, que, en pleine analogie 
avec les autres spectres discutés, la raie 998 K existe en outre de la raie 989 K de la 
pyridine pure et que leurs intensités relatives dépendent de la concentration de la 
pyridine (Mierzecki 1957). Un phénoméne analogue a été observé aussi par Chior- 
boli (1957). = 

En comparant les fréquences du tableau 5 on peut remarquer, que 1) La différence 
entre la fréquence de la raie de la pyridine pure 989 K et la fréquence de la raie nou- 
velle de la solution ne dépend pas de la masse de la molécule du solvent et que 2) 
la raie nouvelle apparait quand la molécule du solvent contient l'hydrogéne actif, 
ainsi que lorsque la molécule du solvent peut former des liaisons de coordination 
(Tramer 1958). 

Dans les spectres des solutions de la pyridine dans l'éther éthylique de l'acide 
acétique on n'observe aucuns changements dans cette région (fig. 2). La molécule 
de l’éther ne contient pas d’hydrogene actif et ne peut former de liaisons de coordi- 
nation, ۱ ; 

’apparition de raies nouvelles dans cette région a été expliquée comme la con- 
séquence de la formation de complexes entre les molécules de la pyridine et celles du 
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7 
250 
e plaque N°55 
o plaque N°58 
200 


x plaque N°51 


150 


100 


100 80 60 40 20 0 


% mol CsHsN 
Fig. 4c 


Fig. 4. Intensités intégrales des raies a) 989 K, b) 1003 K, c) 1029 K en fonction de la concentration, 
dans les spectres des solutions de pyridine dans l'acide acétique. 


solvent. Parmi les différents facteurs pouvant donner comme effet des changements 
de fréquence on peut considérer les suivants: 1) Alourdissement de l'oscillateur (dans 
ce cas la fréquence doit étre abaissée suivant la masse), 2) le couplage des oscillations 
de la molécule avec celles de la liaison intermoléculaire (théoriquement on attend l'aug- 
mentation ainsi que la diminution de la fréquence et on n'observe que la diminution), 
3) les changements structurels. Le fait que la raie correspondant aux pulsations du 
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100 


50 


0 
100 80 60 Zo 20 , 
% mol CHN 


Fig. 5. Intensités intégrales corrigées de la raie 1003 K en fonction de la concentration. 


Tableau V 


Fréquences des raies nouvelles dans la région de 1000 K. dans les spectres des solutions équimolaires 
۱ de la pyridine 


Solvent Formule Fréquence Auteur 

Alcool éthylique C,H,OH 9951 Chiorboli (1957) 

Eau H,0 996 . Krishnamurti (1931) 
Bayard (1943) 
Chiorboli (1957) 

Pyrole C,H,N 998 Chiorboli, Manaresi (1954) 

Chlorphenol OH : C,H, Cl 1000 Chiorboli (1957) 

acide n-butyrique C,H,COOH 1003 Mierzecki (1957) 

acide acétique CH,COOH 1003 ۱ Mierzecki (1955a) 

acide formique HCOOH 1008 Mierzecki (1957) 

Chlorure d'antimoine SbCl, 1000? | Tramer (1958) $ 

dioxyd du soufre SO, 1006 Tramer (1956) 

nitrate d’argent AgNO, 1006? Tramer (1958) 

chlorure d'arséne AsQ], 10154 Tramer (1958) 

iod I, 1005 (infrarouge) Rice, Walsh 
(1952) 


Mulliken et coll. (1953) 


1 Cette raie a été observée pour la premiere fois par Hatem, Valladas-Dubois, Volkinger (1949) 


D'après nos mesures les résultats de l'auteur italien sont plus précis. 


2 Concerne la solution — 2596 mole du solvent. 

3 Concerne la solution concentrée du complexe dans la pyridine (1 
quence dépend de la concentration. 

4 Concerne la solution — 6096 
a un complexe équimolaire. On a observé auss 
spond a un complexe d'une autre composition. 


5% mole du solvent). La fré- 


mole du solvent. D'aprés Tramer (1958) cette raie correspond 
i une raie de fréquence 1000 K qui probablement corre- 
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noyau dans les molécules de pyridine qui forment les complexes a une fréquence plus 
élevée que celle de la pyridine pure doit être alors expliqué par les changements 
structurels des molécules de pyridine. Il n'y a aucune raison pour admettre, que les 
angles du noyau sont changés dans le complexe, car la liaison intermoléculaire est 


V=601K 


v:1573 K 


8, 


Fig. 6. Oscillations du noyau pyridinique correspondant aux raies discutées. 


formée par l'azote qui a une paire d'éléctrons libres. L'hybrydisation de cet atome sp? 
est alors la méme dans la molécule libre et dans le complexe. 

La fréquence plus élevée des raies discutées a été expliquée comme l'effet des 
changements de distribution de la densité éléctronique dans la molécule de la pyri- 
dine (Mierzecki 1957). Cette distribution influence les constantes de force des liai- 
sons et ceteris paribus — les fréquences. Ces changements sont provoqués par le 
caractére éléctrophile de la molécule du solvent et ne dépendent pas de la présence 
de l'atome de l'hydrogéne dans la molécule. 

En comparant les résultats des mesures de l'intensité des raies 989 et 1029 K 
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(fig. 4a et 4c) on observe que l'intensité de la première raie diminue plus vite que celle 
de la seconde. La raie 989 K disparait dans les spectres des solutions ddr. moins 
de 25% mole de pyridine, tandis qu'on observe la raie 1029 méme dans les spectres 
des solutions trés diluées de la pyridine. Sur cette différence entre les raies Ju dou- 
blet pyridinique on a basé la supposition concernant l'origine de la raie 1003 K. On 
a supposé, que cette raie est une raie du complexe et qu’elle forme une eine 
de la raie 989 K (Mierzecki 1957). On ne voit aucune raison d'interpréter cette ed 


-1029 K 
-989 K 


1 i 66.5 % mol 
$ CsHN 
787 % mol 


Fig. 7. Courbes photométriques pour la région de 1000 K des spectres des solutions aqueuses dela pyridine. 


comme la modification de la raie 1029 K. À notre avis la molécule de pyridine pure 
effectue des oscillation d'autres fréquences, que la molécule de pyridine formant 
un complexe. Évidemment une seule molécule n'effectue pas en méme temps les 
oscillations normales d'une molécule de la pyridine pure et celles d'une molécule en 
complexe (Schwab et Glatzer 1957). 

L’intensité de la raie 1003 K passe par un maximum pour la solution équimo- 
laire (fig. 5) (Mierzecki 1956). Ce résultat a été confirmé par Schwab et Glatzer (1957). 


Dans les spectres des solutions équimolaires de 2,6-lutidine on observe dans : 
1007 K environ en outre de la raie de la 


` 


cette région une raie nouvelle de fréquence 
substance pure de fréquence 994 K. 

b) La raie 1020 K. 

Dans la région discutée o 
Dans les trayaux précédants on a consi 
de cette raie à savoir: 

On a supposé que 1) c'est une raie 1 
des solutions (Fénéant 1954), 2) que son appar 


n observe encore une raie nouvelle de fréquence 1020 K. 
déré les trois hypothèses concernant l'origine 


acide acétique renforcée dans les spectres 
ition peut étre liée avec l'existence des 
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ions C,H,NH*, ou bien 3) qu'un autre complexe donne origine à cette raie (Mierzecki 


1955a, 1956). 

Afin d'expliquer l'origine de cette raie on a étudié les spectres des solutions spé- 
cialment deshydratées de la pyridine dans l'acide formique à basse concentration de 
pyridine. Les résultats des mesures de fréquence sont introduits dans le tableau 4. 
Dans les spectres des solutions contenant 7,9 et 209/, mole de pyridine on n'observe 
dans cette région qu'une seule raie nouvelle de fréquence 1010 K. La raie 1022K 
n'a été observée que dans le spectre de la solution composée de 359/, mole de pyridine. 
On a soupçonné que malgré nos précautions cette solution contenait des traces d'eau. 
Ce soupgon a été confirmé aprés avoir comparé les spectres des solutions auxquelles 


1011 K 
1023 K 
1029 K 


n 
۹ 


Fig. 8. Infiuence de l'eau ajoutée à la solution de la pyridine dans l'acide formique contenant 35% mole 
de pyridine, a) traces d'eau, b) 20% mole d'eau, c) 50% mole d'eau. 


on a ajouté de différentes quantités d'eau. Comme le montrent les variations des 
courbes photométriques sur la figure 8, le rapport de l'intensité de la raie de fré- 
quence 1020 K à celle de la raie de fréquence 1010 K augmente à mesure qu'on 
ajoute de l'eau à cette solution. 


La raie de fréquence 1020 K n'est pas alors une raie d'acide, car dans le spectre 
de l'acide formique il n'existe aucune raie dans la région de 1000 K. Elle présente une 
modification d'une des raies de la pyridine. L'apparition de cette raie dépend de la 
présence de l'eau dans les solution. Cette raie pouvait étre observée aussi dans les 
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spectres des solutions plus diluées de la pyridine dans l'acide formique, quand on 
a ajouté une quantité suffisante d'eau. 


c) La région de 1600 K. 

Dans la région de 1600 K on observe les variations suivantes: Dans le spectre de 
la pyridine pure on observe deux raies assez fortes de fréquence 1573 et 1580 K et 
une raie faible de fréquence 1596 K. La raie 1573 K correspond à l'oscillation 8b de 
la symétrie B, (Kline, Turkevich 1944, Corrins et coll. 1953), la raie 1580 K — a l'oscil- 
lation 8a de la symétrie A,, la raie 1596 K — à l'oscillation qui est la somme des 
oscillations 6a et 1 de la symétrie A, des fréquences 601 et 989 K. Cette dernière raie 
est renforcée par la résonance de Fermi avec la raie correspondant à l’oscillation 8a. 
On peut alors supposer, que sans résonance de Fermi la fréquence de la raie 1580 K 
serait un peu plus élévée. Dans les spectres des solutions on observe deux raies d'une 
intensité égale de fréquences 1575 et 1596 K. Cet effet peut étre expliqué de la ma- 
nière suivante (Mierzecki 1957): 

La fréquence de la raie correspondant à l'oscillation 1 est plus élevée dans les 
molécules de la pyridine qui forment des complexes. Les conditions qui permettent 
la résonance de Fermi n'existent plus. La raie de fréquence 1596 observée dans le 
spectre de la pyridine pure et dans les spectres des solutions ne correspond pas dans 
ces deux cas aux mémes oscillations. Dans le spectre de la pyridine elle correspond 
à la somme des oscillations 6a et 1, dans les spectres des solutions — à l'oscillation 
8a d'une fréquence plus élevée. Des observations et explications analogues ont été 
dernièrement publiées par Chiorboli (1957). 

Dans les spectres des solutions de 2,6-lutidine on n’observe pas de variations 
analogües dans cette région, car dans le spectre de la lutidine la résonance de Fermi 
n'existe plus. 

d) la région de 3000 K. 

Dans les spectres des solutions de la pyridine dans les acides saturés on observe 
dans cette région une raie nouvelle de fréquence 3068 K environ. Comme on a montré 
dans le travail précédant (Mierzecki 1955a) cette raie apparait en plus de la raie 3055 K 
de la pyridine. Dans les spectres des solutions plus diluées de la pyridine dans l'acide 
formique on observe encore une autre raie nouvelle de fréquence 3096 K environ. 
Cette raie devient plus forte, quand on ajoute de l'eau à la solution contenant 35% 
mole de pyridine. Elle est observée comme une seule raie dans cette région dans les 
spectres des solutions deshydratées plus diluées avec l'acide formique. Cette raie 
a été observée par Chiorboli (1957) dans les spectres des solutions de la pyridine dans 
les acides acétique et formique à basse concentration de la base. 

Dans les solutions équimolaires de 2,6-lutidine on observe dans cette région 
une raie de fréquence 3067 K environ, un peu plus élevée en comparaison avec la 
raie de la base pure 3061 K. Ze va 

La raie de fréquence 3055 K dans le spectre de la pyridine correspond à l'oscil- 
lation plane symétrique des liaisons C—H. Le ان‎ de cette Se 
les spectres des solutions peut étre expliqué comme l'effet du changement de la distri- 
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bution éléctronique dans toute la molécule de la pyridine, les liaisons C—H com- 
prises. D'aprés Chiorboli (1957) la raie de fréquence 3096 K correspond aux oscilla- 
tions des liaisons C—H dans les molécules ioniques de la pyridine. 


V. Conclusions 


Dans les spectres des solutions des bases pyridiniques dans les acides saturés on 
observe l'apparition de raies nouvelles en plus des raies des bases pures. En particu- 
lier on observe des raies nouvelles dans la région de 1000 K et dans la région de 
3000 K. 

En se basant sur les mesures d'intensité de la raie nouvelle de 1003 K et des raies de 
la pyridine 989 et 1029 K dans les spectres des solutions de cette base dans l'acide acé- 
tique en fonction de la concentration on peut conclure, que 1) la raie 1003 K présente 
une modification de la raie 989 K et que 2) le complexe qui se forme dans les solu- 
tions étudiées et qui est responsable de l'apparition de cette raie, est composé d'un 
nombre égal de molécules de la base et de l'acide et probablement contient une molé- 
cule de chaque substance (Mierzecki 1956, 1957). Lakshmanan (1954) et Chiorboli 
(1957) qui n'ont pas mesuré l'intensité de cette raie donnent d'autres formules pour 
le complexe. Il semble, que leurs résultats ne donnent pas de preuves suffisantes 
à leurs conclusion. ; 

Dans les spectres des solutions étudiées à haute concentration de pyridine on 
observe la raie de fréquence 1720 K (Mierzecki 1955 b). Une raie analogue a été ob- 
servée aussi par Chiorboli (1957). D'aprés l'auteur italien cette raie doit étre liée avec 
les oscillations de la liaison C — O de la molécule de l'acide prenant part à une liaison 
d’hydrogene. C'est pour cette raison, qu'il suggére, que le complexe se compose d'une 
molécule de pyridine et de deux molécules d'acide. 

Syrkin et Sobezyk (1956) qui ont étudié les moments dipolaires des solutions 
discutées proposent pour le complexe une formule I, d'aprés laquelle une seconde 
liaison d’hydrogène peut se former entre le groupe C = O de l'acide et le groupe C—H 
en position 2 du noyau pyridinique. Un noyau du chelat qui se forme ainsi, doit 
stabiliser cette formule. On pourrait expliquer par cette formule l'apparition de la 
raie 1720 K au lieu de la raie 1745 K caractéristique pour la liaison C — O libre. 
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On peut d’autre part expliquer l’apparition de cette raie par l’influence de change- 
ment de la distribution de la densité éléctronique dans toute la molécule de l'acide. 
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Par le cons effet on peut expliquer aussi, que la raie qui correspond aux oscillations 
de la liaison C—C dans la molécule de l'acide acétique est deplacée a une fréquence 
plus elevée (Mierzecki 1953, 1955b, 1957, Chiorboli 1957). Dans ce cas la formule II 


peut représenter la structure du complexe. 


On a montré, que l'apparition de la raie 1020 K dépend de la présence de l'eau 
dans les solutions. Les solutions dans les spectres desquelles on l'avait observée 
n'étaient pas probablement deshydratées d'une maniére suffisante. Il est bien possible, 


que des complexes triples composés des molécules de la base, de l'acide et de l'eau 
sont à l'origine de cette raie. 


L'auteur tient a exprimer a MM. les Professeurs dr J. Pniewski et dr T. Skalinski 
les remerciements les plus sincéres pour l'intérét qu'ils ont porté à ce travail. 
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ON THE THEORY OF MAGNETIC ANISOTROPY OF FERROMA- 
GNETIC CRYSTALS 


By LEON KOWALEWSKI 
Department of Theoretical Physics, Adam Mickiewicz University, Poznan 
( Received June 1, 1959) 


Using the second quantization method of Holstein-Primakoff and Tiablikov, the 
author has calculated: 

1. the magnetic field value and direction dependence of energy of crystals with no 
spin waves excitation, 2. the magnetic field value and direction dependence of dispersion 
formulae, 3. the temperature and magnetic field value and direction dependence of free 
energy and of crystal magnetic anisotropy constants. These calculations were performed 
for simple, body centered and face centered cubic lattices as well as for simple hexagonal 
and closely packed hexagonal lattices. 


1. Introduction 


Ferromagnetic single crystals show anisotropy of magnetic properties. The free 
energy of ferromagnetic single crystals is a function of the angles between the direc- 
tion of spontaneous magnetization and the crystal axes and the direction of external 
magnetic field. 

In the following calculations we shall neglect the spontaneous magnetostriction 
and all external tensions. The anisotropic part of free energy will be then determined 
by the crystallographical symmetry of lattice. 

This part of energy we will call the magnetocrystalline anisotropic energy. 


Akulov (1929) has indirect defined the so called tensor of ferromagnetic aniso- 


tropy using the symmetry properties of cubic crystals. This tensor permits to find 


the formula for the anisotropic magnetocrystalline energy: 
RENE) 
jf oz Ko + Ko + K, %1 %2 ۰ 
2 2 
w = a? a + ai و۵‎ + 2 pa (1) 


between the vector of spontaneous magnetization 


where a, are cosines of the angles 
and the crystal axes, and Kı, K, the so called anisotropy constants. The formula (1) 
ry considerations (Vonsovski, 


can be obtained also by means of rather simple symmet 
Shur 1948). # 
(59) 
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By means of the tensor of ferromagnetic anisotropy we can treat all the funda- 
mental problems of the theory of technical magnetization for magnetic fields causing 
rotation of the spontaneous magnetization. 

On the basis of classical theories we cannot explain the dependence of the aniso- 
tropy constants upon magnetic field and temperature. 

In the paper quoted above Akulov has shown that the magnetocrystalline aniso- 
tropy is caused by magnetic coupling. 

Strictly considered, the magnetic interactions are relativistic effects. Unfortu- 
nately, the relativistic quantum mechanics of many particles is not yet sufficiently 
developed and therefore we have to introduce the magnetic part of energy on the 
basis of correspondence considerations. 

The magnetic energy is very small in comparison with the exchange energy. 
Nevertheless, the magnetic energy represents a dominant factor for the phenomena 
of technical magnetization. This is caused by the fact that changes of energy of lattice 
stationary states, which correspond to isotropic interactions of atoms, and are caused 
by the magnetic interaction energy, are large in comparison with the differences of 
energy for the nonperturbed lattice states. 

Many authors have calculated the influence of magnetic interactions by means 
of the usual perturbation theory. This procedure is incorrect because, as mentioned 
above, the changes of energy levels caused by magnetic perturbation are not small 
in comparison with the differences of unperturbed levels. 

To overcome this difficulty one can introduce a very strong external magnetic 
field, so as to enlarge the differences between the neighbouring energy levels of lat- 
tice. This has the drawback that one can not investigate then the dependence of free 
energy on the magnetic field direction, because in a very strong magnetic field the 
resultant magnetization is directed always along the magnetic field. 

The correct calculation of anisotropic free energy can be performed by the 
methods elaborated by Bogoliubov and Tiablikov (1949) or by Holstein and Prima- 
koff (1940). 

One has to remark (Vonsovski et alii 1955) that the magnetic energy consists of 
a quasiclassical part (Holstein, Primakoff 1940) and an exchange part (Tiablikov 1950). 

In this paper we shall neglect quasiclassical multipole interactions (the long range 
interactions) because, these interactions do not play any larger part in the FES 
of magnetic anisotropy. 

The magnetic exchange energy is caused By not radially symmetric shape of 
electron-clouds at the lattice nodes, In the electrostatic field of ions the electron-clouds 
orientate themselves with regard to the crystallographical axes. The magnetic spin 
moments ,feel" the symmetry of crystallographical lattice owing to their interac- 
tions with orbital magnetic moments. 

For the following calculations we take into account only the dipole-dipole and 
quadrupole-quadrupole magnetic exchange coupling. 

We shall proceed to calculate the anisotropic part of free energy of ferromagnetic 
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single crystal with hexagonal and cubic structure, and we shall obtain the tempera- 
ture dependence of anisotropy constants, as well as their dependence on values and 
direction of magnetic field with regard to crystallographical axes. 


2. The Hamilton Operator 


We assume that the crystal consists of magnetic ions, each with a resultant spin SS. 
For simplicity we shall omit the interactions of magnons with magnons and with 
phonons. 

It is convenient to perform the calculations in the spin space because, from 
the polar theory of metals (Bogoliubov 1949), for spin excitations only, the 
Dirac — Hamiltonian results for each spin value. Fherefore in my paper the 


Hamiltonian is following (Van Vleck 1937, Holstein a. Primakoff 1940): 


SET 2J (Rin) $, 5, x pou (Rin) G F EN Rin) (S, Ry) Rl] + 


+ $2100) S RS, ` Rim)? Ri — 


lm e 


T 


d (2) 


Me 


1 


In the expression above 

J (Rm) stands for the exchange integral er l-th and m-th atoms; 

Rn for the distance between l-th and m-th atoms; 

P(R,) for the coefficient of magnetic pseudodipolar or "anisotropic exchange" 
coupling, whose approximate value is ~ J (R,,,) ° (g — 2)? (Van Vleck 1956); 

Q (R,,) for the coefficient of magnetic pseudoquadrupolar or "anisotropic exchange" 

coupling whose approximate value is ~ 1/4 J (Ry) ۰ (e — 2)*; 

for the "Landé-factor, supposed isotropic; 


for the spin operator of the Lions in h/2 units; 


for the external magnetic field; 
for the value of Bohr-magneton; 
for the total number of atoms. 


دم بیط BB‏ 


3. Eigenvalues of the Hamiltonian 


; CC 
It is convenient to introduce instead of spin-operators the Bose operators bj, b; 
(see Appendix I) by means of relations: 


= 0 — 2 3 + ۵۴۰۶ (ii) by + OF BF (hy 


ro = (1-3). | x 
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Here 
Sf are the a components of c-spin-vectors of lenght ورگ‎ by, b, stand for the creation 
v. z. annihilation operators. 

The کی‎ satisfy the relations: 


5 (592 = S, (Lee OBEN: (4) 


@ 
Q, vectors are defined by formulae: 


Qr = Ps a= eft (S; + Sf) + eml (S, — Si), 


oj = È- S ع]‎ (S, + SP) + ee (S, SH) 
OF = Sr? [S? — (SM) 
et = (SF + iS]) - SD? + (NIT, 6) 


and they satisfy the following relations: 


E OF = Sr [St — (SD, OF OF =. — (5271. 
OF" OF = Sj! [S] — e (6) 


When we substitute (3) into (2) we get in the free magnons approximation 


ques و‎ E (7) 


where 


By = — 4 32) (Ru) Cae i میم و‎ (Ru) SER (S, RRA 


am” یم‎ Q (Rig) Gy Rig)? En Big)? Big 9 3 Q (Ra) + 
lam 
e ES Nes 
AEE p> Q (Rin) Rime (Sp ۰ R)* — u D H- Sı, (8) 


EK 
g — S? — 28 + 1, 


represents the energy of fundamental lattice state in the first approximation. 

The zero-th approximation represents the classical Hamiltonian. 
ur He are small operators of the first v. z. second order of magnitude with regard to 
operators b, and br. 


We shall calculate 5, from the condition of minimum of energy E, with additional 
conditions (4). 


On the Theory of Magnetic Anisotropy of Ferromagnetic Crystals 


We get the following equations: 
N A N 3 NNS 
= 2 2J (Rig) رد‎ + Paes (Ry) [Si — 3 (Sj Rin) Rin Rime + 
N > -> 2 
O) p Q (Rim) (5j: Rin) Ry, Ringe + 
=1 


N 
2 > > m - > > > 
T S? e y. Q (Rim) (S, 1 Ry? (S,, : Ry) Rin Rie 1x uH E m 5 9 


m 
1-1 


where پر‎ are the Lagrange factors. 


Taking into account the formula (9) we can transform #, to the form: 


> 


^ N + 
= > Vin Sm n Tu 
m=1 
where 
eem ar a* p* 


On the basis of conditions (6) we have Hy =0 
The full Hamiltonian has then the form: 


= Ey + Y AQ E, Y B (Rig) 9 di + D C (Ram) bi dm + 
Im lm lm 
+ E D (Rim) bb + 3] E (Rin) Pr bf + را‎ 6 (Rin) bibr» 
Jm lm Jm 


where 
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(9) 


(10) 


(11) 


(12) 


A (Ry) = — AJ (Rig) S + 4۳ (Ry) S — 6۳ (Rn) (QF * Rig) O, Ri) Ri? + 


2 = c: Q (Rim) ) Sr Rim) (S, Ry) (9j - Ri») (2, . Ruy) Bei , 


B (Ry) = ss = (— 352 + 38 — 1) Q (Rin) (GR (Sun Am! Rin + 25? Q (Ry) + 


T E ESS 1۱ یر(‎ 6 Rm)? Ri? + 2 (5 — 1) 0 (Rim) (S, Rn)? Rz, 


CU = = BIRD (A: Rin) (Q, Rin) Rz? + 
Ir a o Q (Rin) (S Ry) (S, $ Ry) (2 E Ri») (2%, E Ry) mon 
DR) = IP (Rim) O Ron) (Om Ry) Ri + 


at E o Q (Rim) (S, ; Ry) (S, ; Rip) (9, x Ry) (Q, & Rim) Rio 
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E (Ry) = O (Bry) (QF ۰ Rig)® (Sop Ri Rint + 
OR) (Qj ۰ R,)* (Q, Rin) (95, * Rig) Rind 
G (Rip) = Q (Ry) (Qj: Ri)? (Sp Ram)? Rif + 
+ Q (Rip) (95 Ron)? (D, - Ron) (95, Rum) Ri. (13) 


We can transform the Hamiltonian above to diagonal form. For this purpose 
we shall pass from space-representation to momentum-representation, or, expressing 
this in another way, from creation (or annihilation) operators for reversed spins to 
creation (or annihilation) operators for spin waves. 

In order to do this, we put 


b = N> JJ exp (— i K,- R) 5,, 
A 


bji NAS) exp ( Kı R) 5j, (14) 
S 
where K! = ze Ki SE Kes s Ay; Ag, Ag assume any integral values 
1 2 


between — 4 G, and 4 G, — 1, — 4 G, and $ وم $ — ,1 — وم‎ and $ Gg — 1, respe- 
ctively. Gi, وم‎ and G, are the lengths of the specimen in the 1, 2, 3 directions in 
units of lattice constant. 

Neglecting the influence of boundaries of the specimen, we obtain thus the follo- 
wing form of Hamiltonian: 


a 


ê = E, + J U, j bı + DV} + ۶ P (15) 
À 
where 


U, = YMA(R)exp EK, R) + B(R)], R,— R,— Rp» 
Lead 


Vi = DIC (RD exp (i K,- R) + E (Ry), 
Pr = DG exp (i K,- R;) GRAN (16) 


stands for the summation over all the next neighbours of any atoms.‏ زر 
h‏ 


If the magnetic interactions Vf and V7 were neglected, the Hamiltonian above 
“would be diagonal. The Hamiltonian (15) is not diagonal, because operators b, b. 
and b, b*, do not commute with 5} b,. 
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In order to diagonalize the full Hamiltonian we use the methods of Holstein and 


Primakoff or Tiablikov (Bogoliubov 1949, p. 206) putting: 


ak ak 


À 1 2 A A 
bı = gg exp (ip) [a (êa + êd + (8 لس‎ 


l 4 A A * * 
Yep (ipa) [h (ĉa — 623) + klärt En]. (17) 


We define this transformation for the positive halfspace of K , only, i. e. we put 
K a = — K,. 


We confine our attention to the case of |K,| <1, i. e. to the long spin waves. 
Then 


ide E Set Ce ۳ U, — )-۳ Û 
مار‎ TEINTE AREA 


By use of the formulae (17) we get the following form of the Hamiltonian: 
D 


where کم‎ 


pp = 2 (E, — Uy), (19) 


Ax A 
and c}, c, are Bose operators. 
" ^*^ 
The eigenvalues of the operators c, c, are n; = Os + 


4. The Fundamental Energy Levels of the Crystal at 0*K 


We assume that the crystal in the lowest energy state is magnetized to satura- 


tion, i. e. the S, vectors do not depend upon the ۰ 
The cubic lattices: 


E, Eat Kyo — Np HS. 1 (20) 


We get the second anisotropy constant K, if we take into accounte the hexapole in- 


teraction. 
a. The simple cubic lattice: 


+ = > 
R,=R rp = "Th 
des > am مت‎ 
پر‎ = E ep E ey cei eee 


a — lattice constant; " E e, stand for the reciprocally orthogonal unit vectors. 


# 
5 
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0 OOO ne رآ‎ Eee ee 
There are six nearest neighbours for each lattice node. In the zero-th approximation 
(E®) i. e. for the classical Hamilton function we get: 
Ko = — 6NJS? + NQS*, K; = — 2NQS*, 
In the first approximation (Ej) which corresponds to a state with no spin waves: 
K, = — 6NJS? + NQS? (S? + 2), 
KE, = — 2NQS c. 


We get from the conditions (9) for energy minimum the following values for 
spin-components: 


I. 

2 2 
for St== 0: “SY = ipis CENE 
for S90: "= ME, Ses 
for 9 20: bx ids benc d 

i. e. crystallographic directions of the type [110]. 
II. 


ege Sa 0 Se E. 
for S*— $*—0: PS, 
for See mu 0: LS = E E 
i. e. crystallographic directions of the type [100]. 
III. 


og eu? AN 


geb و شیب و‎ us 


i. e. crystallographic directions of the type [111]. 


In the zero-th approximation we get therefrom the results shown in the table 
below: 


Table I. 


Direction 


[100] — 6NJS? + NQs* maximum value of EQ) | minimum value of EO 
1110] — 6NJS? + 1NQS* intermediate value intermediate value 

= — of EO of EO 
[111]- - — 6NJS? + 1 NQS* minimum value of EQ maximum value of EO 


e‏ و 
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In the first approximation we get the values shown in the table IT. 


Table II. 

Direction Ep S+1,0>0 S+1, Q<0 
[100] — 6NJS? + NOS + 2NQS? maximum value of Ej | minimum value of E 
[110] — 6NJS? + 1 NQS? (S?+2S+ | intermediate value intermediate value 

+3) of Ey of Ey 
[111] — 6NJS? + 1 NQS? (S?+4S-+ | minimum value of EA | maximum value of EQ 
+4) 


For S — 1 the energies for the three main directions are identical. 

It can be seen that the direction of easy magnetization depends on the sign of the 
coefficient of quadrupole interaction Q. 

b. For the body centered cubic lattice (e. g. iron, S = 2): 


_ 3 
pr 


Ry QE ê + Èy + 69), Ra 


E 3 E E E 
> (esc: e, + êz). 


There are eight nearest neighbours for each lattice node. In the zero-th approxima- 
tion (EQ) we have: 
K, = — 8NJS? + 4/9 NQS*, K, = 16/9 NOS), 
In the first approximation we get: 
Ko = — 8NJS? + 1/9 NQS? (45? + 16S + 5), 
Ki = 16/9 NQS? c. 


In the zero-th approximation we get the values shown in the table III. 


Table III. 


3 0 
. [100] — 8NJS? + 4 NQS* minimum value of EO maximum value of ES ) 
[110] — 8NJS? + $ NQS* intermediate value intermediate value 
of Ej of Ef 
A o. 
[111] — 8NJS? + 28 ۶ maximum value of Eo minimum value of ze 
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The values for the first approximation are shown in the table IV. 


Table IV. 


Direction Ey SLOS 0 
[100] — 8NJS? + 1 NQS? (4S? + 16S + 5) minimum value maximum value 
of Ey of Ey 
[110] — 8NJS? + 1 NQS? (8S? + 85 + 9) intermediate value | intermediate value 
of Ey of Ey 
[111] — 8NJS? + 4 NQS? (28S? + 16S + 31) maximum value minimum value 
of Ey of Ej 


For S — 1 the energies for the above three directions are identical. 

The direction of easy magnetization depends on the sign of the coefficient of the 
quadrupole interaction Q. 

c. For the face centered cubic lattice (e. g. nickel, S ~ 1): 


R= > Hate, hA. F(+4+%), 


There are twelve nearest neighbours for each lattice node. In the zero-th approxima- 
tion we have: 


> = — 12 NJS? + ۷۵94 K, = NQS*. 
In the first approximation we get: 
K; = — 12NJS? + NQS? (8? + 2S + 2), 
K; = NQS? c. 


For the zero-th approximation we get the values shown in the table V. 


Table V. 


Direction Q<0 


[100] — 12NJS? + NQS* minimum value of EI) maximum value of EO) 
[110]. — 12NJ$? + 3 NQS* intermediate value intermediate value 

| of EQ of EQ 
[111] —2NJS*--$NQS* | maximum value of E® | minimum value of EQ 


ems 
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For the first approximation we have the values shown in the table VI. 


Table VI. 


Direction E Ej | S+1,0>0 S+1,Q0<0 


[100] — 12NJS? + NOS? (S? + 2S + 2) minimum value maximum value 
of E] of E] 

[110] — 12NJS? + 1 NQS? (55? + 6S + 9) intermediate value | intermediate value 
of Es of Ej 

[111] —12NJS? + 1 NQS? (4S? + 4S + 7) maximum value minimum value 
od Ey of Eq 


For S — 1 the energies for the above three directions are identical. 

The direction of easy magnetization depends on the sign of the coefficient of the 
quadrupole interaction Q. 

The hexagonal lattices: 


E= a + b (893 + e (594 — Nu H- S 
= aga je sin? à + K: sint à — Nu F- S, 
Ky = ab S CS Ky = (b He’ St) St, eye n ATE 
$ — is the angle between the direction of spontaneous magnetizatiori and the 
hexagonal axis. 
a. For the simple hexagonal lattice: 


T = 2 s IR DI 
Ry = Rath = daîn, Th =k 6x, cE €x; ۳ 


+ 


€x: ee e, stand for the reciprocally orthogonal unit vectors; €, is the unit vector 


directed along the hexagonal axis. There are eight nearest neighbours for each lattice 


node. In the zero-th approximation (EQ) we get: 
Ko = — 8 NJS° + NPS? + NQS*, 
K, = — 3/2 NPS? — 2NQS* < 0. 
K, = 17/8 NQS* > 0. 
In the first approximation (E, we have the values: 
Ki = —8NJS? + NPS? + NQS* (S! + 5+ 2), 


Ki = — 3/2 NPS* — OS (252 — 3S + 1), 


Ki, = 17/8 NQS? c. 


I 
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If only the six nearest neighbours in the basis plane are taken into account, the 
easiest magnetization direction is perpendicular to the c-axis. 

Moreover, when we take into account also these two neighbours, which are on 
the c-axis below and above the basis plane, we find that the easiest magnetization 
direction will still be perpendicular to the hexagonal c-axis, but the minimum of 
crystal energy for this direction is not as deep as in the former case. 

b. For the closely packed hexagonal lattice: 


8 
R, = Rit ip == GA, y Fe 


There are twelve nearest neighbours for each lattice node. For the zero-th approxi- 
mation (EQ) we get: 


K, = — 12NJS? + 4/3 NQSS, 
— 2/3054 > 0, K, = 7/12NOQS*>0 


A 
I 


For the first approximation (E) we have: 

Ki, = — 12NJS* + 2/3 NOS? (28? + 25 + 5), 
K, = — 2/3 NQS? c, 

Ki = 7/12 NQS? o 


When we take into account all the twelve nearest neighbours for the closely 
packed hexagonal lattice, we find, that in comparison with the case of the simple 
hexagonal lattice the direction of the easiest magnetization leans over toward the 
hexagonal axis. Whether the easiest magnetization direction lies along the c-axis or 
deviates from this axis is determined by: S — the spin value, the influence of more 
distant neighbours and the hexapole interactions. 


5. The Influence of Magnetic Field on the Spontaneous Magnetization 
Direction at 0°K 


When no external magnetic field exists, the spontaneous magnetization is directed 
along one of the easiest magnetization directions. When the external magnetic field 


H forms the angle ©, with the direction of easy magnetization [hy ky Lol, the spon- 
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. taneous magnetization will lie in the plane determined by [ho ko lo] and H, formin 
an angle © with the direction ] ko Lol- | i 


fh, lo] 


Fig. 1. Directions of easy magnetization, spontaneous magnetization and external magnetic field. 


O increases with ۰ 
From the condition for minimum of energy E, we obtain the dependence of © 


upon the value of H: 


dE, 


dO 


Eo FE Des = Jale I * cos (Oo RA 0) 


at 


we get thus: 


d E anis. 
dO ; : (21) 


ies Ve sin (Oo pee 0) 


The magnetization in H direction is: 
I = I." cos (0, — O). (22) 


From (21) and (22) we obtain the 1 = I (H) dependence. 

For example, when we take [ho ko Lol = [100] and H is directed along [110], 
the direction of spontaneous magnetization in a single crystal of iron will coincide 
with that of magnetic field when the value of field reaches 345 Oe. 

The magnetic field influence on the rotation of spontaneous magnetization was 
phenomenologically investigated by Bozorth (1936), (1951), Schlechtweg (1936), 
Tarasov and Bitter (1937). 


6. The Dispersion Formulae 


The value of E, (19) or approximately of U, represents the energy of spin wave 
with wave vector EK In the calculations below we consider the nearest neighbours 


only. 


y 
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a. For the simple cubic lattice we have: 
U, = A + a? B, (K3)? + a? B, (K3)* + a? B, (K2 (23) 
Can 
A = a ao + E H-S, 


B, = b+ b ai + b ai, 
B, = b + b, a + b, ۰ 
B, = b + b, a2 + b, «à, 


a, = 8Q5* (1 — S), 


ag = 8QSn, 

n =5S?—7S + 2, 

b = (4J + P)S, 

b, = — (3P + 400) S, 
by = ۰ 


In the present paper the exchange integral J and the pseudodipolar coefficient P 
have values four times less, the pseudoquadrupolar coefficient Q is sixteen times 
less and the spin number S two times larger than the corresponding expressions used 
in the paper of Pal (1955). It follows, that for the magnetic field directed along the 
spontaneous magnetization, the above mentioned formula is identical with the Pal’s 
formula. 


b. For the body centered cubic lattice (e. g. iron S ~ 2) we get: 


U, = A + a? B (Kj)? + 2a? C (K K3) + 2a? D ) K3) + 2a? E (KY K%), (24) 


A =a + موه‎ + 7 3 


I 


AIS — 2/9 050 + 16/9 QSooc, 
— 4PS 02 &, 

— APS a a, 

SAPS qi. 

a, = 6405 (— S? + 4/3 5 — 1/3), 
dy = 64Q8n. 


SMS E, Cu 
۳ of 


/ 
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c. For the face centered cubic lattice (e. g. nickel, S ~ 1) we obtain: 
U, = A + a? B, (K$) + a? B, (K3)? + a? B, (K3)? + 2a? C K3 K3 + 
+ 2a? DK} K3 + 2a? EK} Kj, (25) 
A =a, +a,0 + RES, 
D = by + b, o1 + b, at + bs oi, 
B, = b, + b, 2 + by و‎ + 5, 0, 


B, = by + b, ۵ + ba ۵ + bo 
C = — (3/2 P + Qo) S a, a, 
D = — (3/2 P+ Qo) S 0 «3, 
E = — (3/2 P + Qo) S «a ag, 
a, = 4QS (S? -- 17S +2), 

,40 — = وه 

b, = 1/4 (16JS + PS — 2QSo), 
b, = — 1/4 (3PS + 4QSo), 

b, = — 1/2 QSo, 

QSo.‏ = و 


d. For the simple hexagonal lattice we have: 


U, = A + a? B (K3)? + a? C (K3)? + a? D (K3)? + 2a? EK; KY, (20) 


2 ER > 
A =a, + aso + ag aê + HS, 


B = b, + ba وه‎ + by وه‎ + ba ai + bs a1: 
C = + c « + cg وه‎ + c4 at + cs ۰ 
D = d, + d, a3 + d; ag, 

E = e + وه‎ a? + eg a3, 

a, = 3PS + 1/2 QS (— 9S? — S + 6), 
a, = — 9PS + QS (39S? — 47S + 14), 
a, = 1/2 QS (— 855? + 1195 — 34), 

b, = 6JS + 3/8 PS — 3/8 QSo, 

b, = 9/8 PS — 3/4 QSo, 

b, = 9/8 QSo, 

b, = — 9/4 PS — 3050, 


EI? 
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b; = 3050, 

c, = 6JS — 15/8 PS — 9/8 QSo, 
c = 27/8 PS — 3/4 QSo, 

Cs = 21/8 Ose; 

c, = 9/4 PS + 30So, 


1 
ce = — 30050, LE | E | 
` d, = AJS + PS, fe à RT 


d, = — 3PS — 40So, 

4QSo,‏ = و 

e, = — 1/2 OP + 600) S a, a. 
وه‎ = 3/2 0So a, €, 

QoS a, a.‏ 9/2 = وه 


Ebr the case of axial symmetry we get: 


= A + a? B (K? + aè C (KY? + aè DG, : (27 


2 
A 9 


Sc 


B =b, + baag + baj, me FE 
Ci a re T a 
D SEN D 5d + daa + dy ai, pn Ge? 3 5 ER d 
EN TU de 
where dy, ag, وه‎ have the values as above. Dt. o — | CG 


de 
FE 
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e. For the closely packed hexagonal lattice we get: 
U, = A + a? B (Kj)? + a° C (K) + a? D(KJ? + 2a? EK; K3 + 
+ 2a? FK1 Ki + 2a? GK} Kj, (28) 


De 
A = a, + وه‎ 05 + as وه‎ + As 


B = b, + b, al + b, ag + ba a? + bs at + bg a a3, 
C — ده‎ c a1 + c, و۵‎ + cy وه وه + وه‎ + 2 5 
D = d, + dy o + dg وه‎ + d4 %3 + ds ai + de «f o. 
E = e, + و6 + %2 وه‎ o5, 

F = fi fa %3 fat 

G = gi + وق + %3 وق‎ 23> 

a, = — QS (S? +55 — 6), 


da = — 205m, 

as = 5/3QSn, 

b, = — BIS — 3/2 PS — (8 + 7/12) Q3o, 
b, = — 1/2 PS + 58/3 Qo, 


b, = — (10 + 5/9) QSo, 
b, = SPS + 64/3 QSo, 


b, = — 12 QSo, 

b, = — (19 + 1/6) 050, 

e, = — BJS — 15/2 PS + 1/4 QSo, 
c, = — 1/2 PS + 5/6 QSo, 

cy = — 23/36 QSo, 

c, = 5PS — 3/2 QSo, 

cg = — 73/18 QSo, 

cg = 9/2 QSo, 


d, = — 8JS — 6PS + 19/27 QSo, 
d, = — 6PS — 2/27 QSo, 

d, = 31/27 QSo, 

d, = 8/27 0, 

dy 
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d, = — 8/27 QSo, 

e, = 1/2 (P + 1/3 Qo) Sa, ae, 
وه‎ = 1/9 Qo Sa, ae, 

ey = — 7/6 Qa Sa, xs, 


fi = (4P + 7/3 Qo) Su 2; 
fs = — 2/3 Qo Sa zs 
fa = — 29/9 Qo Sa, Gs, 


gy = 4/3 (P + Qo) Say a5, 
وق‎ = 32/27 Qo Suez «3, 
gg = 4/27 Qo Say ۰ 


For the case of axial symmetry we obtain: 


U, = A + a? B (K$)? + a? C (K3)? + a? D (Kj)? + 2a? FK; Kj, (29) 
A =a, + as od + مه‎ ad + A 5 
B = b, + b, a3 bs a’, 

C. m tait cg ai, 

D = d, + d, وه‎ + ds ۰ 

F —fi- fao, 

a, = س‎ QS (S? + 5S — 6), 

ge = — 2 QSn, 

as = 5/3 QSn, 

b, = — 8JS + 7/2 PS + 3/4 QSo, 
b, = — 11/2 PS + 22 QSo, 

by = — (3 + 7/18) QSo, 

€ = — 8JS — 15/2 PS + 1/4 QSo, 
وه‎ = — 1/2 PS + 5/6 QSo, 

c4 = — 23/36 QSo, 


d, = — 8JS — 15/2 PS + 1/4 QSo, 
d, = — 6PS — 10/27 QSa, 
dy = 39/27 QSo, 
۹ f = (4P + 7/3 Qo) Sa, «3, 
fa = — 29/9 Qo Sa as. ` 
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Coefficients of K}, K3, Kī consist each of an isotropic and an anisotropic part. The 


anisotropic part, which is of magnetic pseudodipole and pseudoquadrupole origin, 
amounts to about two percent of the isotropic one. 


7. Free Energy 
Free energy of crystal is given by formula (Appendix II): 


pe kT J, mfi- )هه‎ (30) 
À 


Going over from the sum to the integral (Appendix III) we get 


+00 
N 1 
F= Et | nf — ep (7) (A + a? B, (Kz)? + 


+ a? B,(K3)? + a? B, (Kj)? + 2a? CK; Kj + 2a? DK; Kj + 
2a? EK} “|| 06 ۰ 010 ۰ 0 و‎ )31( 
where N, is the number of lattice sites in one simple cell. We transform the above 
quadratic form into the canonic one. 
a? [B, (K3? + B, (K3)? + B, (KD? + 2616 KX + 2016 K + 2EKT Ki] 
= a? [y (KE)? + Xe D2 Fs (D3 (32) 
where y,, X» and xg are the real roots of equation: 
B, — و‎ G D 
G By — x: E ves 0, 
D, TE ss! Va 


The coefficients of the transformation, which leads to the canonic form of (32) are 
given by the following equations: : 


B. me Xo C; D Bix 
C, By uc E E Box — 0, 
D, E, B, aT Ak Pax 
thus 
Kx Pi Bor وم‎ [RÀ Kj 
KY | = | Bio Boo Bs: Ke BI 


ES Bis Bos Bsa! MS Ki 
# 
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The inverse transformation will have the following form 


Kj Pu Ba Pis (Ki 
K3 = Baı Kë Bas K , 
A: Pai fsa Bas ۶ 

When we substitute: 
X1 À x së € 5 d = Xs ۰ z => 
ya Kj =u, a K} =v, a ET Ki =o, 
we get 
5 NB 
F= E, (kT. 
o + (AT) N, Gn? 
where 


-L (y), 33) 
۰ à? Vos Xe Xa = 


+00 


L(y = [ma eut: 


(u? + v? + w?)]} du ۰ dv do, with y = ITO 
Introducing spherical coordinates o, £, Ë and putting oi = ۶ we get: 


L (y) = 27 | In [1 — exp (— v — 9] Vt dt 
0 


Performing the integration (Appendix IV) we get 


3| e exp (— ny) 3 

rw 7 (5) وت رز‎ - -2 E ۰ ۸ (y), (34) 
We calculate approximatively the A(y) sum using Euler’s method (Bieberbach 1930 
Valenta 1955): 


TESTS 


SEEN ۳ (4) dx + + [g(eo) + g 0([ + 
E Le (es) — 2 (I+ fev (x) ۰ W (x) ۰ dx, 

| ۱ 1 

where 


x NNI ند‎ 
gE (=x ‘exp (— yx), B = 
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After some calculations we obtain: 


A (y) = (11/8 — 5/4 y) exp (— v) + 4/3 o Va [1 — % (V y)] + R (y), 
where 
5 Yv 
۵9 - ۰ | exp (—n2) dp, 
0 


is the Gauss integral (Jahnke, Emde 1938). 
We proceed now to the estimate of R(y): 


Hi 1 = f f 
RE [ie EI de PNE مد > مره "وال‎ IE 
1 n= 1 1 


= 10-2 (y? 4 IS: — 3 : 5/2 ۲ y? 2 Lg 3 ; 5/2 x 7/2 WI I 5/2 : 7/2 v 9/2 Lay) 
= 107? (y? + 5y + 35/4) exp ( v). 


where 


m 9 
ENTE Im= — emie) "E 
2 2 m —2 2 


1 


For y <0.1 the rest R(y) represents maximum 6.5 percent of A(y) value. 
Using the formula (33) we shall now calculate as an example the free energy of 
simple cubic lattice (for y <0,1): 


F = 4 + B (kT) + Č (ET) + D (ET) + 
+ [E ET)” + 6 (e) + DR cos (S, H) + Ko, 
Ka 4 ET + GT + ET + 
+ [e (KT + E (KT) ] H cos (S, H), 


A = —6NJS* + NQS? (S? + 2), 
B — 11/8 ۰۵, 

C ۸۰6۰۵۰ مره‎ 

D ۸۰6۰ E: a,n 

E 2۰۰6۰ 

6 —2u-:a: 4۰۵۰ E, 
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Nc A 
a = — NOUS! 
B =11/84-¢, 
y =A- yt an 
à — LAC USE Ge: 
= = 2۸۰۸۰2۰۵ 
E =u: A TS, 
"us NEEL 
(27)? a3’ 


$m eus (a 5 v 


E = 5/3 ya [1 — 9 (Vw)l + 31/16. 
y = 1/24 (4 Va [1 — ۵ (Vy)] — 63}, 
t = 1/4 (4JS)-"* - S? (OP? + 64JQn). 


Appendix I 


Let us introduce two orthogonal systems of axes: 

l. a system x; y, z which is rigidly bound with the crystallographical structure of 
lattice; 

2. a system x’, y’, z' whose z'-axis is directed along the direction of sponta- 
neous magnetization. 

We can express the spin operators in 5, representation in h/2 units by Bose 


operators b,, b;, 
Sf = YS, [f Gi) 6, + b; f (0J, 
P Aay’ S A N * ^ 
= px Ua) 6; — b f (ĉa), 
so 35) 


‘The transformations (35) are known from the paper of Holstein and Primakoff (1940). 


The ری‎ S", S” operators are matrices, which transform as simple polar vectors. 
The spin operators in the system of coordinates x, y, z are: 


Sa cos (x, x") SN + cos (x, y) SY + cos (XE Se, 
S? = cos (y, x’) S* + cos (y, y’) SY + cos (y, z) 5%, 
S* = cos (z, x’) S* + cos (z, y) SY + cos (z, z") S. (36) 


On the Theory of Magnetic Anisotropy of Ferromagnetic Crystals 81 


When the A axis lies in the xy plane, i. e. when X (z, y) = x/2 (see Fig. 2), 
we obtain for S*, S", and S* the formula (3). 


Fig. 2. 


In the A units these formulae take the form: 
Oa a ny 1 a ^ e il ek Ax X 
S= SF [1 — چ‎ + = Qf f (nj) bı + 5 Qr di f (fu), 
Si 2 2 


When SF tends to 5j, S? = 0 and S7 tends to zero, or when X (y, y’) = 9, x (x, x’) > 0 
and À (z,z’) > 0 we get transformation formulae of Holstein and Primakoff. 


5f = = Lf (tu) bi + bi Fi), 


yn KG LS (ê) b — b f (AD), 


Az Ak 4 A pie Ty 
Sf = Sı — by ` br, f= 1-5 


Tiablikov (1956a a. b) has also obtained the transformation (3) using a method of 


approximate second quantization. 
It can be seen that the transformation 
Tiablikov are the special cases of (36) formulae. 


formulae of Holstein — Primakoff and 
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Appendix II 
Free energy is given by: 
F= —kTln Z 


where Z is the sum of states. 


Z = 5» exp | = d (nra x Em) | 
1 


{74}=0 


— exp E (E; + [هه‎ 7] exp [= em Ey 2] 


470 


1 oo 
ONT | 77 He + هه‎ D 


FEM AE LIT 3 In D - اجه‎ Sall 


4= -0 


Appendix III 


We can replace the summations >} which extend from — $ G; to 4 G; — 1 
(14343) 


by summations from — oo to + oo, because the short spin waves induce only negli- 
gible corrections to the value of free energy. The function which is summed, depends 
on the scalar products, which are invariants for transformation to another coordinates 
system. 

Instead of summing over the indices 24, Ag, Ag, which describe any translation 
lattice with cubic or hexagonal arrangement, we shall sum over A,, Ay A» which 
correspond to a simple cubic lattice. In regard of the very high values of G,, G,, G,, 
the summations over À,, Ay, A, can be replaced by integrations: 


=. GrGyGe [os 32 
Dr = is jo. 


We calculate G, G, G, from the following condition: 
GG, G: = 
> EI DS 


where the integration is extended over the first Brillouin zone. For cubic lattices 


cw 
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G, G, G, gives the number of simple cubic cells = N/N,. N stands for the overall 
number of lattice sites, and N, for the number of sites lattice in one simple cubic cell. 
For hexagonal lattices there is a fractional value of N,. 


Appendix IV 


To calculate the integral: 


L' e | ple 
(y) = J ep en dt. 


we consider the more general integral: 


co co 


| = a= fry exp [— n (y + 01] dt 
H 0 


n=1 


oo 
= [rare ay) ap (9 dt, 
n=1 9 


Putting nt = و‎ we get therefrom 


A exp (— r =~ exp (— ny) 
year for exp (— ٩ dq Fil Y, SES, 
0 n=1 


n=1 


Hence 


This is the (34) formula. m 
I should like to thank Professor Dr. S. Szezeniowski for helpful advice, detailed 


discussion and direction of this work. 
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ELECTRICAL RELAXATION IN NI—ZN FERRITE 


By L. 7 
Department of Ferromagnetics, Institute of Physics of the Polish Academy of Science, Poznan 
( Received June 15, 1959) 


The electric and magnetic properties of iron deficient Ni-Zn ferrite were investi- 
gated to obtain an explanation of the mechanism of dielectric relaxation. The relaxation 
observed was found to be related to electron diffusion, and. not to inhomogeneous 
structure. 


Introduction 


The electric properties of ferrites are as yet insufficiently well known, notwith- 
standing the fact that these are materials of great importance for basic investigations 
and practical application. Their unusually high dielectric constant and its dispersion 
are of especial interest. A number.of papers (Polder 1950, Koops 1951, Fairweather 
& Frost 1953, Volger 1954, Yoffe et al. 1957, Peters & Standley 1958) deal with these 
properties. The first four give an explanation of the dielectric properties of ferrites 
based on the inhomogeneous structure of polycrystalline ceramics. Koops attempted 
to explain the dielectric properties of Ni —Zn ferrite by assuming the presence of two 
layers within the grains of the polycrystalline material. The equivalent RC system 
yielded frequency dependent resultant capacity, losses and resistance in accordance 
with the curves of these quantities as obtained by measurement. The theory was gene- 
ralized recently (Parker & Smith, 1958). It would seem, however, that agreement with 
the theoretically predicted results has been due more to the choice of definitions 
(Hippel 1954, p. 88), whereas in reality the phenomena observed were caused by rela- 
xation effects detected at about the same time and related to electron diffusion (Kan- 
ichi-Kamiyoshi 1951). This is suggested by some doubts as to whether the assump- 
tions.on which Koops has based his computations are justified. According to his theory 
very thin boundary layers*) would have to surround grain cores of greater conducti- 
vity. However, it is well-known that the rise in conductivity of ferrites is related to the 
appearence of Fe** ions (Wijn 1953, Domenicali 1950, Calhoun 1954) as a result of 


* Some authors place the thickness of the boundary layer at about 1 À (Volger 1954). 
(85) 
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partial reduction in the process of heating which should be greatest at the surface of 
the grains. It would seem that re-oxidation during cooling should be impeded by the 
fact that the Fett ions are chemically bound. If such is the case in reality, additional 
difficulties arise for Koops theory, on account of his assumption, that the dielectric 
constant of the boundary layers has the same value as that of the grain cores. This 
assumption is contradicted by results cited by the same author and showing that both 
the conductivity and dielectric constant depend strongly on the Fett ions contents. 
In this context it is interesting to mention also results by Yoffe, Khvostyenko and 
Zonn, who proved i. a. that, in CO—Zn ferrite, the properties of polycrystalline 
ceramics resemble those of monocrystals. Finally, Peters and Standley carried out 
measurements of the dielectric properties of Mg—Mn ferrites that did not corrobo- 
rate Koops' results. 

The present investigation is the first of a series aimed at explaining the mechanism 
of dielectric relaxation in ferrites. 


Results 


ZnO „_ NiO,, Fe,0 


amples. 
0,7 0,5 comp 


30,3 

The samples were prepared using standard techniques, by heating a compressed 
mixture of oxides at 1200°C and letting it cool gradually with access of air. The grain 
size was about lu. X-ray analysis, carried out by Mr Przedmojski of the Chair of 
General Physics B of the Warsaw Polytechnical School, showed the samples to contain 
ZnO in addition to material of spinel structure. The lattice constant amounted to 
8.40 A. The method of preparation was found to influence the properties critically. 
The dielectric and magnetic properties were measured by the resonance method 
using an espécially adapted Ferisol French-made Q-meter. A general diagram of the 
measuring device is given in Fig. 1. The dielectric measurements were carried out with 
the condenser c, especially adapted to use throughout the range of — 60° to + 300°C. 


Fig. 1. Simplified diagram of the measuring device 1. — Magnetic stabilizer, 2. — Power supply, 
3. — 50 kc — 70 Mc oscillator, 4. — Vacuum-tube voltmeter. 
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Fig. 2. Temperature dependence of dielectric losses at different frequencies. 
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Coating with electrodes of any kind (graphite, silver) was found to present no advan- 
tages in the case considered. Checking by dielectric measurements and by measuring 
the conductivity (using a Wheatstone bridge of high precision) proved that heating 
the sample to about 300°C produces considerable changes in its properties after some 
time. This could be avoided by using plate electrodes pressed onto the sample. 
Although in the latter case considerable difficulties arise if the exact absolute values 
of the conductivity, dielectric constant and losses are to be measured, carrying out the 
measurements with all possible precautions and care yields the thermal dependence 
of the foregoing properties with good reproducibility, which is especially important. 
Moreover, plate electrodes held pressed to the sample reduce total losses by lowering 
the contribution from conduction, thus enhancing the degree of precision of the 
dielectric measurements. Measurements were made for different temperatures at 
constant values of the frequency. The permeability was measured on a toroidal sample 
with a single copper winding. In this case, the dielectric losses were neglected. An 
ultrathermostat was used for investigating both the electric and magnetic properties; 
however, the samples did not contact the silicone oil used in thermostat, as it was 
found that even traces of oil in the sample bring about considerable changes in its 
properties. 

Typical curves of the losses are shown in Fig. 2, wherein characteristic maxima, 
shifting towards higher temperatures as the measuring field frequency rises, are seen. 
From a number of such curves, Fig. 3 was obtained, showing the relation between the 
frequency and the temperature at which maximum losses occur. Curve 2 was obtained 
with a sample additionally heated at 360°C during 5 hours. The anomaly about 
240°C is related to the Curie point, as proved by measurements of the initial perme- 
ability. The results are shown in Fig. 4. Maximum magnetic losses occur at the tempera- 
ture at which the permeability drops to 0.8 of its maximum value. Measurements at 


0 100 200 300 T°C 
Fig. 4. Temperature dependence of the initial permeability at 10 Mc. 


Electrical Ralaxation in Ni—Zn Ferrite 89 


other frequencies yield identical results. For comparison, the measurements of the 
electric conductivity shown in Fig. 5 were carried out. The latter reveal a disconti- 
nuity at 155°%. At the Curie temperature there is no such discontinuity. Curve 2 
was obtained after the sample had been heated at 360°C for 5 hours. The discontinuity 
is displaced towards higher temperatures. Curve 3 shows measurements on a sample 
that had been additionally polished down to about 50% of its original thickness. 


108- 
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2 25 1000 
Ty 
Fig. 5. Temperature dependency of the electric conductivity. 


This proves that if the results are to reproducible, uncontrolled polishing cannot be 
employed. A rise.in the activation energy shows the deficit in oxygen to be related to 
the sample surface. 

A measurement of the thermoelectric power proved the ferrite to possess electron 


conductivity. F: 
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Discussion 


The foregoing results cannot be explained by Koops'theory. Although it is diffi- 
cult to give a correct interpretation in the case under consideration, it can be shown 
that even in so complicated a case important information may be deduced. 

In principle, the dielectric relaxation observed in polycrystalline ferrites might be 
interpreted as arising from inhomogeneous structure (Wagner, 1914). Eq. (40) of the 
paper cited yields the following characteristic relaxation time for a system consisting 
of a non-conducting medium of dielectric constant ¢ wherein small spheres of a mate- 
rial possessing the same dielectric constant e and the conductivity x are suspended: 


3 € 6 
Tas dx (1) 
Although Wagner's assumptions are but very roughly approximated by Koops'- 
model of ferrite ceramics, by assuming x = x, exp (> St it is neverthelees possible 


to obtain the grossly correct temperature dependence of the relaxation time: 


E 
ER Lä 2) 


However, a comparison of curves 3 and 5 points to the circumstance that the 
temperature variations of the point at which maximum losses occur are not related direc- 
tly to the conductivity, since not all points of discontinuity coincide. This is a hint that 
relaxation is not related to inhomogenity of the material. In addition to inhomo- 
genity, the relaxation observed may reside in a displacement of the ions or electrons. 
To consider this possibility, it should be stated that Figs. 3 and 5 may serve for 
determining the activation energies for relaxation and conductivity. They suggest 
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Octahedral Tetrahedral 


Zn ++ Z ++ 
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Fig. 6a. Ion distribution over all possible lattice positions. 
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` It is noteworthy that E; = & Ez ~ €»: 
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that the relaxation must be related to electron transitions between ions of different 
valency, since the conductivity is of electron type and the activation energies differ 
but very little. Fig. 6a illustrates ion distribution over all possible lattice points. The 
dashed line denotes antiparallel spin array (negative exchange integral), the con- 
tinuous line-parallel array of the spins (positive exchange integral) (Néel 1948). 
Figs. 3 and 5 show that the exchange energy of a pair of ions in the relaxation pro- 
cess must be positive, whereas that of a pair of ions participating in conduction must 
vanish. The change in activation energy below the Curie point is probably related 
to the presence of structure-sensitive conduction (Koch and Wagner 1937; de Boer 
and van Geel 1935; Nijboer 1939). The displacement of this anomaly towards hig- 
her temperatures on heating the sample points to the participation of Fe** ions in 


Lon ER 
۳ are HE 
Fe Fe 
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Fig. 6b. Proposed scheme of electron jumps which might explain observed phenomena. Single 
line — relaxation, double line — conductivity. 


both relaxation and conduction. On heating, their number diminishes because of 


the diffusion of oxygen into the sample (Smolenskiy 1952), and clearly the con- 


tribution from defects and impurities begins to play a part. 


The. foregoing makes it plausible that the most probable scheme of electron 
trasitions is the one illustrated in Fig. 6b. It should be stressed that the possibility 
of conduction related to electron transitions between ions of iron in the tetrahedral 
positions has been dismissed, as modern super-exchange theory does not exclude 
the possibility of J44 Z 0 (Weisz, 1952). 

The curves of Figs. 3 and 5 can be described by the following eqs.: 


E E 
o= A, exp EI + وم‎ (a) 


f= B, exp LA + Byexp LA 8 


The constants A,, and Dua may be computed by extrapolating to infinite tem- 
perature. Thus, e. g. Bı and B, should be of the orders 1011 and 1013, respectively. 


Another interesting point is that the mechanism of magnetic relaxation is not 


identical with that of dielectric relaxation, a fact resulting from Fig. 4. 
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It should be concluded that, at least in the case under consideration, Koops’- 
theory alone is unable to account for the dielectric properties of ferrites, though 
the existence or absence of boundary layers has not been decided definitely one way 
or the other. A good description of the phenomena is obtained if electron diffusion 
is assumed to be the cause of relaxation. 

Results of investigations on other ferrites will be published in the near future. 

The author wishes to express his thanks to Professor Dr S. Szczeniowski for 
his continual encouragement of the present investigation. 

The author is also indebted to Professor Dr. S. Krupicka for his stimulating 
discussion, and to Mr. J. Baszyñski for his help in carrying out the measurements. 
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GENERAL LAGRANGIAN AND CANONICAL FORMALISM DE- 
SCRIBING THE RELATIVISTIC MOTION OF A FREE ROTATING 
PARTICLE 


By F. HALBWACHS 


Institut Henri Poincaré, Paris 
(Received August 5, 1959) 


The present paper is an endeavour to formulate the problem of the relativistic 
rotating particle in a Langrangian form by availing oneself as far as possible of the main 
results of a previous hydrodynamical theory of the author. 


After several attempts to deal with this subiect (Aymard 1956, Unal and Vigier 
1957; Takabayasi 1958), we have recently set forth, together with P. Hillion and 
J. P. Vigier (1958), a Lagrangian treatment of the relativistic hydrodynamics of 
a fluid endowed with intrinsic rotation, the state of this rotation being described 
‚by local Einstein —Kramers variables, i. e., by a local rotating tetrad of unitary orth- 
ogonal four-vectors. Starting from a well-known variational procedure commonly 
used in classical field theory (Belinfante 1937; Rosenfeld 1940; Winogradski 1956 
and 1959; Halbwachs 1959 — Appendix C), we derived from the usual tensor den- 
sities suitable hydrodynamical quantities fulfilling appropriate equations (Halb- 
wachs, Lochak and Vigier 1955, Takabayasi 1957, Halbwachs 1959 — chap. IV). 
In this procedure one uses a Lagrangian density and its derivatives with respect 
to the relativistic coordinates x". Such a formalism (as well as all the work hitherto 
done on spin fluids), has a very abstract and postulatory character, because it is rather 
difficult to think of a proper or intrinsic rotation in a continuous fluid. It appears 
thus paradoxical that the much simpler and more concrete problem of the relativistic 
rotating particle was not yet treated in a Lagrangian form. Such a treatment is just the 
purpose of the present paper, where we endeavour to carry over the main results 
of our previous hydrodynamical theory into particle treatment, using a scalar Lag- 
rangian and its proper-time derivative along the world line followed by the particle. 


I. Basic vector variables 


We shall start here from the idea that a rotating particle is necessarily one 
with some extent in space, possessing thus an infinite number of degrees of freedom 
which describe all the details of its stucture. In order to express the motion of such 

E (03) 
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a particle by a finite number of equations, a certain number of global quantities 
must be introduced by integration over the volume of the particle, the evolution 
of these global quantities being bound by dynamical equations. The more detailed 
the description, the more numerous global quantities we shall hawe to introduce. 

We shall base all our considerations on the two following assumptions: 

1. All the global quantities to be introduced must be relativistic vectors. 

2. A "center of the particle representing its mean position can be defined. 

The particle will thus be described in its internal structure in a more or less 
detailed manner by a system of vectors localized at its center, in the same way as 
a function may be described by its Taylor expansion, which accounts mor or less accu- 
rately for the form of the function in an extended region by means of successive 
derivatives considered at one point. We shall attempt in these conditions to derive 
the laws of evolution of the vectors from a suitable Lagrangian formalism and then 
to build the corresponding canonical formalism. 

In a recent paper (Halbwachs and Vigier 1959), we hawe shown that a Lagran- 
gian founded on the above assumptions and describing a simple particle on which 
no external field is acting, namely ZG IP, qP) (x, denoting the relativistic coor- 
dinates of the center of the particle, a! — the vector quantities localized at that 
point, and the dot standing for the derivative with respect to the proper time 7 of the 
center), yields always two general laws which can be regarded as conservation laws 
and which are the consequences of the scalar character of the Lagrangian, the prin- 
ciple of stationary integral being taken into account along an arc of the world line 
followed by the center. We have namely 
Conservation of linear momentum: 


Conservation of total angular momentum: 
Mw = E X, — G, x, 


These laws, which hawe been known for a long time (Frenkel 1926, Weyssenhoff 
and Raabe 1947) do not depend on the form of the Lagrangian, nor on the kind 
and the number of the vectors d"! used. The vector G, and the antisymmetrical tensor 
M,, are expressed by means of the Lagrangian by the formulae 


o£ 
Ga = uL 
EE 


5) 9, d 
244 ag? 


where the summations are to be performed both over the lower tensor indices and - 
over the upper index r labelling the vectors. These expressions are very abstract, 
and the form of the Lagrangian shall be regarded as adequate only if it leads to giving 
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them some physical meaning, so that eventually G, and M,, will appear as standing 
for the physical linear and internal angular momenta. 

The main problem to be solved in order to arrive at the Lagrangian is to find 
the expression of the kinetic rotation energy term. We shall therefore introduce two 
kinds of physical global quantities, the first kind expressing the rotation in a kinem- 
atical way (the global angular velocity), the other kind in a dynamical way (the global 
angular momentum). The vectors so introduced will be localized both at the “center 
of matter density“ defined by Bohm and Vigier (1958) in their theory of rotating 
fluid Tnasses. 

In order to define the kinematical rotation, we shall first place ourselves in the 
proper frame of the center of matter density at the instant v. Let us assume that 
the extended particle may be looked upon as constituted of a certain number — 
which may be even infinite — of discrete identical material points, which will desc- 
ribe bare matter, irrespective of the energy (for instance, if the particle is electrically 
charged, these material points will be assumed to carry equal electric charges and 
we have to deal then with the distribution and current of electricity). If we then 
study the spatial distribution of these points and of their instantaneous velocities at the 
instant 7, we are led to define a system of three orthogonal rotating unitary vectors 
po, po, p, in such a way that the whole of the material points stands on an average 
at rest with respect to these vectors, so that the instantaneous rotation of this system 
of vectors at the instant t describes globally on an average the rotation of the matter 
of the particle (Halbwachs 1959, chap. V). 

The projections of these three vectors on three fixed axes in the proper space 
make up a matrix Cj fulfilling evidently the relations 


C, G ES du 


Its time-derivative will stand for the instantaneous rotation state of the moving triad 
in the proper space. 

The relativistic instantaneous motion will be completely described by the super- 
position to the above spatial rotation of a Lorentz transformation from the laboratory 
frame to the proper frame, which includes therefore a relative velocity v. Setting 


we shall have a unitary velocity u, 
; er 2 
Uj = Gi U4 = Gaio n E A, 


Along this unitary velocity, we shall take the unit vector 


H = u,/ic, 


that is 2 
Hl = av, lic, bY? = ۵ 
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(It is an imaginary time-like vector, with a positive magnitude, imaginary space- 
-components and real time-component in the Minkowski frame). 

If we consider the three rotating unit vectors BO, £2, go i in the proper space, 
we see that they are the projections on the fixed axes in the proper space of three 
world vectors bP, bo. b®, whose projections on the vector b are zero, and which 
thus constitute Mokether with the latter an orthogonal tetrad of unitary relativistic 
four-vectors. 

The projections BY, BY, B® are just the matrix elements Ci mentioned above, 
and the relation Cj C = پر‎ merely assures the orthogonality and unitarity of the 
three vectors in the proper space. 

The expression of the four vectors bo» b, b, be with respect to the labora- 
tory frame will be given then by the classical Lorentz transformation (Takabayasi 


1958) 


= any, bY = ۾‎ 
um j ok @ _ (ck a Vi Uj k 
ba P C Cj ? bj Ci RE (1 + a) c j (6) 


. These formulae can be looked upon as defining the ‘tetrad b® in terms of the 
ordinary space velocity v of the center of the particle with respect to the laboratory 
frame and of the matrix C? which defines in a conventional way the global proper- 
-space rotation. To summarize, the four world vectors H) are the so called Einstein- 
-Kramers variables defining in a kinematical way the Get relativistic rotation of 
the particle. They obey the two equivalent sets of orthogonality and unitarity con- 
ditions 

pO X = Ox > p pm = ô", 


In order to define precisely the global relativistic angular velocity, we start from 
the elementary treatment of the nonrelativistic dynamics for rigid rotating bodies. 
It is well known that any rigid rotating body admist at each instant of an axsis ins- 
tantaneously at rest, and the instantaneous rotation state of the body at this instant 
may be described by an axial vector c lying along this axis, the so called rotation 
velocity vector. Two of the Einstein-Kramers vectors, namely B® and f? considered 


in the proper frame, may be chosen in a plane orthogonal to the instantaneous rot- 
ation axis. 


We shall hawe then 


> 


de NE Ee ^ 
emm = e) (D — DO (2) 
e UN E 


that is, if we denote by a dot the derivative with respect to the time (which in the 
proper frame, is just the proper time of the center of matter density) 


= ۰ 1 n 5 . H 
w; = ëm BY Be = Eir By BE = 3 Eijk (By ^d BP BEP) 
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But it occurs frequently (namely in gyroskopic problems) that it is more suitable 


to choose the (8%, B)—plane to be orthogonal to some axis B® fixed in the rigid 
body, such as for instance the inertial axsis of a top undergoing a precessing motion. 


Its velocity has to be superposed to that of the vectors go and B® and the angular 
velocity is then given by 


ex = à eg. (By? BE? + BP BY? + By By”) 


It is just the latter form we hawe to choose in order to describe the average rot- 
ation of a complex (or eventually fluid) system for which the concept of instantaneous 
rotation axis does no longer belong to the existence of a set of material points momen- 
tarily at rest. We thus write in the proper frame 


Eijk (Br BE),‏ $ = ,س( 
where the summation over r = 1, 2, 3 is understood. If we introduce the relativistic‏ 
four-vectors bo (whose time components are zero in the proper frame)‏ 
0 } = 
(o9 = 3 ej, OP OP).‏ 


Finally, we may add to the sum in the right hand side the homologous term CN 
b(®)0 owing to the vanishing of the vector b in the proper frame. We thus have 


(o)? = $e, Gf P), (E — 1, 2, 8, 4). 


The so called axial vector œ; considered in the three dimensional space is known 
to be the dual of an antisymmetrical space tensor, what is brought out in the above 
expression by the appearance of the symbol Se. In order to pass to the relativistic 
covariant form, we have to look just upon this tensor as stanting for the angular 


velocity, that is 
(alt =} (b? H — be HEI, 
which is simply 
(059 = (© MÉI, 
the latter tensor being antisymmetrical owing to the orthogonality conditions. This 


expression may become relativistically covariant by adding (in the proper frame) 


the homologous term 
(op)? = OP BPP. 
As in the proper frame 


1D = bP = 09 =0, P=) 


d 
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our supplementary term (w,4)° reduces to 


0 
(or = ay = Lay - (2), 
that is to the linear acceleration of the center of matter density. 

If (w,,)° and (هره)‎ are regarded as the components in the proper space of a rel- 
ativistic antisymmetrical tensor w,,, then we see that the latter embodies at once 
the rotation of the particle and its acceleration, both properties being brought out 
as separate components only in the proper frame, while they are inseparably bound 
together in any other frame. 

Finally, we shall introduce as the generalized rotation velocity the relativistic 
antisymmetrical tensor 


o, = 19 9. 

In order to build the complete expression of the relativistic rotation kinetic 
energy, we have to introduce yet another set of vector variables expressing the dyna- 
mical aspect of the rotation. We shall consider a priori only the cases where this 
aspect may be described by an antisymmetrical tensor, namely the proper angular 
momentum S,,, which is just the relativistic generalization of the angular momentum 
of elementary dynamics. Such an assumption was put forward by all theorists who 
endowed point-particles with proper rotation (Frenkel 1926, Weyssenhoff and Raabe 
1947). It is also by such a tensor that the angular momentum of the electron may be 
represented, the quantum average value of which is provided by Dirac’s theory. 
The authors who ascribed to the rotating relativistic particle extension in space (Mat- 
hisson 1937, Moller 1949, Bohm and Vigier 1958) have also put into evidence similar 
tensors. 

In particular, in Möller’s theory the particle is described as a droplet of relativistic 
classical fluid, characterized at each point by a symmetrical energy-momentum den- 
sity t,, the divergence of which vanishes if no external field is acting. Möller introduces 
then the proper angular momentum density 


May, = (x, JE X,) la — (x, 7 X,) TE 


which is just the moment of t,, with respect to the so called center of gravity X „that 
is the barycenter of the droplet relative to its energy distribution t,,. By integration 
of m,,4 over the volume of the droplet, one gets the total angular momentum S,,. 
Bohm and Vigier (1958) follow the same way, but they take as center X, the so called 
center of matter density, which is the barycenter of the distribution of matter (or of 
electricity). They get thus a more general theory. 

Whatever the definition chosen for the tensor S,,, it is always possible to split 
it in a covariant manner into two relativistic vectors by projecting S,, and its dual over 


the unitary velocity (Karpmann and Raman 1956, Halbwachs 1957 and 1959 chap. V) 


LL (4) 
U, = ich, - 


ee P m 
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We thus get the spin 


= 1 1 (4) 
Sn = er U > Ewaß Sap De 


and another vector, which we have called the “dullard” (balourd in french) 


l 


Lm = Kon Si uy = i uv Dë 
These definitions yield for the angular momentum the expression 


Sur = V uva Ze by SS i (08 L, Ax Ox H. 
Obviously 


EE 
5,5, = tp b, = 0. 


This splitting is profitable as it puts immediately into evidence the important special 
case where t, = 0, it is, the well known Frenkel relation S,, u, = 0, which is introd- 
uced as a postulate by most of the authors who have studied the relativistic angular 
momentum (Frenkel 1926, Costa de Beauregard 1943, Weyssenhoff and Raabe 1949). 

Of course the angular momentum, although it describes the rotation in a differ- 
ent way than the angular velocity, is not deprived of any relation with the latter, as 
it depends after all on the internal motions of the matter inside the extended particle, 
which also determine the angular velocity. Thus the tensor S,,, although it was intro- 
duced irrespective of the vectors (S is in fact bound with them through the local 
variables which describe the details of the structure of the particle. We shall make no 
detailed assumption about this structure, but we shall be led nevertheless to postulate 
various kinds of connections between S,, and the bo and this will lead us to various 


forms of the Lagrangian. 


II. Lagrangian formalism 


We shall here confine ourselves for the sake of simplicity to the cases where 


the Lagrangian depends in a linear way upon the derivatives H of the Einstein— 
Kramers variables. As the rotation kinetic energy term is $ میرگ‎ Dia: this assumption 
will imply that the angular momentum depends only upon. the b®, but not upon 


‘their derivatives. 


We shall introduce this dependence by successive steps and in a more and ee 
strict way. As a first step let the expression of مرگ‎ depend only upon the vector b; 


or u,, then, as was shown above 


i E 1 t 
Su, = E Env ap Ue Sp + CH (Uu by — Un Ne Sy: 
ZP ` 
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This expression is just a linear one, and it introduces besides the p two vector 
variables s, and t, characterizing dynamically the rotation. We shall thus assume the 
proper rotation term in the Lagrangian to be 


: Ge 1 TORIO) | 
چم‎ Sue Our = —. E Euvaf Ua Sg + y (Uu t, — u, 7] b® of 
E : l j(r) p (r) 
= cun Mr == t, — ut E 
2 E Euvap Ua Sp + E (Up H A b, b, + 
1 IET 1 à; 


(the label r runs from 1 to 3). 

We note that, in view of the form of the terms in brackets, if the vectors s, and t, 
are split into components respectively orthogonal and parallel to the velocity, the 
latter will vanish because of the antisymmetry. Thus, the time components do not 
appear in the Lagrangian and do not interfere in the motion. We may assume there- 


fore the relations 


to be verified without any restrictions for the motion, and thus we need not introduce 
any corresponding Lagrange conditions. If these relations and the antisymmetry of 
Egg are taken into account, the term with u, u, becomes merely 


f TE 
272 ty uy = — 9c dh 
The other terms to be put into the total Lagrangian are the classical Lagrangian 
} mou, U, and a Lagrange term which assures the orthogonality and unitarity 


conditions, namely 
3 Uy U 
Au RR = [pap toe — 0) 


A, being an essentially symmetrical Lagrange coefficient. We thus get finally the 
following Lagrangian 


Le E Mg Un Un — - B Suva Ua Se + = (Up ty — Uy Al E b? — 
reen E (wor er). ۱ 
ER uy u 9y £5 ۱ 
It is easy to derive the corresponding canonical linear momentum 
Gu = = = Mo Uu — = Soap OS. BY? sy — 


EE Re A ba Ar 
— 9, ty (by b; — by bi) — acl 20 ie 
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In order to derive the canonical angular momentum, we put 


oL re 1 ۲ 1 1 

2 Zum PELA Ug Sp + 17 (u, ۸ n) [b= ES x Bak, 
9L Si 1 

3۳ mite 


and derive therefrom 


oD OL OL eL 
Mr bo — سس اون‎ = — ty 
Or ETS) 27 ^ "aL ER 


= i e ` bO pn (r) pin) l 
= 64 À v =S Dy b) + gc Un t, — Uy ty) 


Ji A 
= Pile + D (Sua u, U, — S, U, u,) + 9c (u, t, — wu, t). 


The last line follows from the last but one by taking into account that bp p = 
= 6, + u,u,/c?. As Spa Ua = ct,, we get finally 


Therefore, the internal angular momentum which is provided by our general 
formalism is identical with the one we introduced in our Lagrangian in virtue of phy- 
sical considerations, so that the proposed Lagrangian is suitable for expressing the 
dynamics of the rotating particle under consideration (with S,, u, 0). 

Nevertheless, the two conservation equations which are, as we know, automati- 


cally fulfilled 


Gi = Ope Se = Cait, Cu 


are not sufficient for the complete determination of the evolution of the thirteen 
variables u,,, Cora a and not even of the twelve variables p DES if we take into 
account the various orthogonality and unitarity relations. (Note. that the second 
set contains one variable less than the first one, because the motion is evidently 
gauge invariant with respect to any rotation of the bP, b in the plane orthogonal 
to the B®, bW). 

Nevertheless, as the Lagrangian formalism provides just as mamy Euler equa- 
tions as the Lagrangian admits of variables, the motion is necessarily determined 
throughout. It is easy to find the necessary additional equations by taking into account 


the Euler equations relative to s, and £,, i. e. 
dL LN 


dO IC NUR 
ALT pe i NÉI uy = 0 
a 


1 i 
B A bP ug = = ES €uva 
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That is 
AR 
Wy U, = 0, 


y v 


T oT ae : P‏ و 
eA pO BO) — BM po GER 0‏ 
Em lT E‏ 


As we have 

u, (69 H — 0۵ bM) = u, 
this yields 

u, (BE HI — 29 HEI = 0, 
that is 


Dy by = 0. 


In other words, the tensor w,, and its dual are orthogonal to the same vector u,. 
This is only possible if the angular velocity w,, is intrinsically zero. So the proposed 
Lagrangian stands for a well defined particular case of the general class of motions 
governed by the system of the ten conservation equations (system of Bohm— Vigier). 
In that case, the purely kinematical (or electrical) matter of the particle, which differs 
from its energy, travels as a whole in uniform motion without rotation and accele- 
ration (if no external field is present). 

Nevertheless such a particle has a non-vanishing proper angular momentum. 
This is not inconsistent, as the angular momentum describes the motion of the energy, 
and in particular the variations of the internal stresses, so that the energy distribution 
may vary with time without any change of matter distribution, as would do for 
instance a rotating heat flow inside a fixed solid, inducing a kinetic momentum 
without any rotation of matter. 

If the Bohm — Vigier motion is made more precise by imposing some restriction 
upon the dynamical variables or by binding the latter more strictly to the kinematical 
ones then, on the contrary, one will get some kinds of angular velocities less strictly 
determined. ; ۱ 

Thus we may assume (Halbwachs 1957, Bohm and Vigier 1958) t, and s, to be 
collinear; this can be expressed by a Lagrange condition 7, Egg Sa tg Or in a more 
simple way by lessening the number of the variables and putting 


S, =O, COS A, t,=o,sin A. 


The Euler equations relative to o, and A provide, after some calculation, 
^ 
(w,, cos À + w, sin À) u, = 0. 


> ۱ 
©, sin A —w,,cos A) U, 0, = 6, u, cos A. 
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We may suppose further (Haibwachs 1959) s, and t, to be mutually orthogonal and to 
have the same magnitude sọ, by introducing the two Langrange terms 


AS and BGS, LL) 


which provide after some calculations 


A AN 

Our Uy = “3 Ua Sp (Wag tu + Mag Su)» 
0 
il A 

Our Uy = — Ua Sp (Wap Su — Wag tu). 


0 


A more interesting assumption consists in imposing the Frenkel — Weyssenhoff 
relation 5,, u, = 0 which implies that, in the proper frame of the particle (or, more 
precisely, of its center of matter density), the only non-vanishing components of the 
angular momentum are the space components. We then get for the thirteen variables 
Uys Gur Su the thirteen equations 2 


Gs 0, Suv = G,u, — G,u S0 — O0 


2 uv Uy 


which form the so called Weyssenhoff system. This system, which was completely 
integrated by Weyssenhoff, yields for the thirteen quantities under consideration 
a well defined law of motion. Nevertheless, it is easy to see that this system is not 
sufficient for describing all what the Lagrangian tells us about the motion, as it does 
not inform us about the angular velocity. ; 

The Frenkel — Weyssenhoff condition can be introduced by adding the Lagrange 
term A, Bun u, . But it is obviously more simple to remove merely in the above Lagran- 
gian the terms containing t, , because t, = مرگ‎ u,/c is now zero throughout. The Euler 
equation for t, will then disappear and only the first equation ر91/95‎ = 0 will remain, 


that is 


^ 
ر0‎ Uy = ۰ 


If we place ourselves in the proper frame, where u? = 0 the latter relation reduces 
to Oo, = 0, that is 


Thus only the time components of the relativistic angular velocity do not wanish in 
the proper frame or, in other words, the matter does accelerate but it does not rotate. 
"We note this interesting complementarity between the properties of the angular 
momentum and those of the angular velocity (Halbwachs. Hillion and Vigier 1958). 
The preceding assumption which ascribes to the angular momentum the expres- 

sion 
Sun B — Ewap Ua SB = — 8 S S8 
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binds only the dynamical variables to one of the vectors of the kinematical tetrad, 
namely AT. We get another kind of motion if we introduce the supplementary condi- 
tion that the vector s, standing for the spin lies steadily along another of the vectors 
beo, let us say for instance b(®. This result can be obtained either by a Lagrange 
condition 
Tu» Nuvap HI Sp 
or, in a more simple way, by setting s, = Ag b where hy = Vs, s, is the magnitude 
of the spin, which is variable in principle. 
We get then for the proper rotation term 


iho 3) 1 
“pig BO AP BY, 


E E 


ae 
The terms with & = 3 and € = 4 are zero because of the antisymmetry, and we are 


left with 


h = 
GE eue DE BY (BY BM + BY BM), 


But we have also between the b® the following relations 
(4) p3) DA) — 2 
Sra Un. DPTO es DOS 


(4) Hä p — (1) 
Carafe Dp O و‎ 


so that we get the expression 


h . ۰ 
x Cr MP — MP P) 


or, because of the orthogonality, 
( (2) 


It follows that the vector DE does no longer appear in the rotation Lagrangian. 
In order to avoid its redundant appearance in the orthogonality condition it is con- 
venient to replace the form used as yet by the equivalent form 


zen (b® im = dëm. 


where we assume & and 7 to take only the values 1, 2, 4. 
Finally, the latter condition being explicitely written down, the Lagrangian be- 
comes 


/ 


1 
L= 2 Mg Uy Uy — hy ou iD en 


a By? — 1) + 429 (BS? — 1) — AM ES tal | 


+ 2412 بو‎ b2 Ti qus Ou uj + "E a - hu, 
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We know already the two conservation relations 


Ban Ee ee 
with 
eL 


OU 


G 


Duy = hy deo pon Ka te U2) 


(the latter expression being tantamount to ho = £,,55 سر‎ bu): We shall have also, 
because of the form of S,,, the Frenkel-Weyssenhoff relation 5,, u, = 0. The whole 
motion, described by means of the dynamical variables, will be the Weyssenhoff 
motion, but the two variables HI), b® which remain purely kinematical ones, do no 
longer fulfil the relations of the preceding case. 

The three variables u,, 6 Be yield the three Euler equations 


A44 A14 724 
Im - 2 د(‎ E DON 2 — bi = Gus (1) 
: (2) (2) a) ALS 
— hig MP + 23:2 D 4-24 +2 — u, = 0, (2) 
TE 167 yeh a SECO, AS 
hig BO + hy BOD + 2212 MP 4 240 02 + 2—— us = 0. (3) 


If we contract (1) by u,, we get 
jM 


Doo c ee RU d 
Gu uy = — 0 273 ۰ 


On the other hand we note that in all equations obtained by derivation, we meet the 
difference m, — 24**/c? and never m, or 45^ separately. It is thus suitable to put 
them together in one term. In order to keep my as the proper mass and to assure 
the unitarity of u,, we shall replace both terms } Mo Up Uy and 44t (u, u, + c?) by 
3 mo (u, uy + ch, mo being regarded as a Lagrange multiplyer which yields a cor- 


responding Euler equation. We thus get 


Aie 124 
š — ain. (1) SE 9 ED pa 1 
so that 

Gele mo C2 


On the other hand, if we contract (2) by b® we see that A1? = 0. Hence, the 

contraction of (3) by b® shows that ho — 0. The magnitude of the spin is constant. 

The value of 411 is derived by contracting (2) by oy and that of 422 by contracting 

GJ by o: 

| 2411 = hy و‎ 62 b = — 2 T, 
4n 949 que a NT: 
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The value of 211 or 422 stands fot the kinetic rotation energy T. If we take into account 
these results, the three Euler equations may be rewritten in the form: 


14 
Mo Uy + 2 = be SE ace TE Gu (1) 
: 14 
me A D TOI ۲ u, = 0, (2^) 
f 24 
ho bo) OT pO = u, = 0. (3) 


We can determine the expression of the two remaining multipliers by contraction 
of the two last equations by u„ which provides 


and these expressions introduced into (1’) give 
ho هن‎ 2) __ (D5 هار‎ 1 
Gu = Mo uy + zz (bà ub, — bj Uab; )= mo u, ۰1 ورگ‎ 


This is a well known relation which can be derived directly from Weyssenhoff’s laws 
of motion. Finally. by contracting (1’) respectively by b® and 6%, we get 


EEE 14 
Gb = لا‎ du = 2 Se 
(2) ho هر‎ AM 
Gu bp = — a bu à li, = 27—. 


We can finish the determination of the laws governing the kinematical variables 
by borrowing new results from Weyssenhoff’s theory. Let us first recall (Halbwachs 


1956) that the spin s, = 0 b is a constant vector, as s, = 0, so that we have b3 =(. 
On the other hand u, is known to fulfil the following differential equation (Mathis- 
son 1937, Weyssenhoff and Raabe pr 

Mo Up SE =m Ip =O) 


If S, is replaced by A, (bP AH — bP b®) and if that equation is respectively 
contracted by HU) and b® we ‘get the system 


Be ER 
Tp tba H oe ti, b® = 0, 


۱ Je: 
mo و‎ d ü, bP = 0. 
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MT MMMM 


waren determines throughout the monu of b® and b because u, is fully deter- 
mined by the theory of Weyssenhoff. It is easy to see that these equations are gauge- 
۰ " ۰ D 5 ۰ B ۰ 1 2 ۰ 

invariant, nos is, they re valid iod and b are respectively replaced by H? 
cos « + b^ sina and — b;? sin æ + bt” cosa, no matter what value is given to ۰ 


III. Canonical formalism 


We shall now use the latter Lagrangian — as the most simple example of a linear 
Lagrangian — in order to build up a relativistic Hamiltonian formalism for the rota- 
ting particle. It is well known that there are two kinds of relativistic Hamiltonians. 
Let us begin with the scalar Hamiltonian, which will yield covariant canonical equa- 
tions. Let us first rewrite the Lagrangian 


= 


L = > mo (Uy u, + c°) — ho KD £2 + 9 ms DDD + 


a 


14 EE 
+ 241 (D BD — 1) + 232 (9D A — 1) وب‎ = bP u, + 2 T Ur wr 


i 


We have obviously as canonical momentum conjugated to x, 


n au pe 
== = Mg Un -+ 2 v og = 2 TER DA 


Ou, 


and as canonical momentum conjugated to AC) 


9L 
(2) 7 d (1) 
x em D Ex is Di. 


Thus 8( happens to appear already in the Lagrangian. For simplicity of notation 
we shall omit from now on some upper indices by putting 


BP = bus UP = — FE Sw = bu Bo = br Bo 


As regards the other variables, as their derivatives do not appear, they provide merely 
the customary Euler equations and their conjugated momenta are permanently zero, 
so that they bring no contribution to the Hamiltonian. The Euler equations would 
have to appear as the corresponding canonical equations. 

Moreover, to simplify writing, we shall put 


124 14 422 Me ae d 
fcc im), etis 
ic o ol SIDA ho 


Thus the Hamiltonian, derived by a Legendre transformation whose generating 


function is the Lagrangian, will take the form 


H= Gin, +B, b, — L(Gys bir Bur mo u PT, TRN) 
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with the canonical equations 


3H aH _; | 
8G, Uus 9B, Hu? | 
3H : aH | 
RE E SS aS 


aH 9 OH _ OH OR) 
$m. . E ye a 


By expressing L in terms of the canonical variables, we get without difficulty 
H = G,u, 4- B,b, - (6, — 2۳ b, - 2B,) u, — 3 ۵ ۵ EE 
ST E h2) — 2 u,b, — 2Lu, p, = 
= pmr FOL, —Ab, B, — Libur —1B,u, — T (b, 5, — 1) — 
— T’ (f, B, — (۰ 


That is, after replacing u, by its expression 
1 
H= pes (C 2l b; 2L Ee mz c2] — 


MOO E TPS ht) — Ab, Pp- 


The canonical equations relative to G, and x, provide 


oH 1 4 H 
86, jams (Gu = 21 b, — 2l B,) Up 


G, = 0. (4) 


D 


This is the expression for the linear momentum and the first Weyssenhoff equation. 
Taking into account the first canonical equation (4), the six remaining canonical 
equations provide the four already known results 


= Fh — =0 (that is „u, = — E 
EH  — 9 ty = 0, Se = — By = 0 

He. man 
En E 
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being the orthogonality and unitarity relations, and 


OH < H ^ ۳ 

Qui r^ d Bac Age bs 

oH : 

à. = — 2 u, — 2T b, — AB, = — By. 
u 


The expressions for the subsidiary variables و وا‎ T, T", A are obtained by con- 
traction of the two latter equations respectively by u,, b,, and f 


4» the orthogonality 
and unitarity relations being taken into account, 


WES 


Aes 


, 1 d 
i} P= ag Êu un 


5 dE ik n 
on ur Nena 
۸ = — bp 6, = 8, B, = 0. 


By substituting these expressions into the complete equations, we get two fol- 
lowing equations 


bg = ببس‎ unt + bf 
139 7 
Bu = SS Ex De Uy * Uy =F p» b, D 


By translation into the tetrad symbolism, these equations become 
(2) LG 4 LO p(D) AO) 
jg — (ED BO) BO + (HP ae A 
1 ; 1 2 2) ° 
iD — de at -D + ÈP BP) «OP 


; : ۵ A 1 ۰ 3 7 (1) 
Thus the derivative b® reduces to its projections on b® and u,, and the derivative bu 
to its projections on b® and w,. The components of both derivatives parallel to by 
A 3 
vanish and thus we are led once more to the result that the spin s, = 0 bo of the free 


Weyssenhoff particle is constant in time. The corresponding equation a = 0 is 
equivalent to the second Weyssenhoff equation. | 
Finally it is useful to write down the Hamiltonian equations in canonical form by 
elimination of the subsidiary variables J, l’, T, T' and À. Equation (4) contracted 
respectively by b, and D, provides 
1 ANN) 
P D Gabe l= 578 ۰ 
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Thus the canonical equations become 


Gan : 
CG, = D Ca Gide bj CEU, UMS 
0 


G : 
G, b, CE Ca by: bu — SN = Bän — bu by), 


G, D, Cr : 
Cobe (e, ann, a). 
0 


This Hamiltonian formalism is very suitable for putting into evidence the cons- 
tants of the motion with the aid of classical Poisson brackets. It is easy to verify that 
we have, besides the constancy of G,, 


[G, u 


p Um SE 


expressing the constancy of the proper mass, and 


[Emap Ur ba Bo > H] = 0 


expressing the constancy of the spin vector (and of course of its magnitude). For 
these calculation u, has evidently to be expressed in terms of the canonical variables. 

We shall finally construct the non-covariant Hamiltonian # leading to non- 
covariant canonical equations expressed in terms of the derivatives with respect to 
the time ۸ = x,/ic. It is well known that # will then stand for the particle energy with 
respect to the considered frame of reference. 

We see easily that if we substitute the above canonical equations in the expres- 
sion of the scalar Hamiltonian H, we get H = — mech, Then, if we write the latter 
equation as a relation between the sixteen canonical variables, we may consider one 
. ofthese variables, namely G,, as being a function of the fifteen remaining ones. We put 
therefore 


eG, X Era 9 D.) 


and get 
9X de o dw o dX — 9x du 28 db, 
9x, di Cd pig a copo MP dE 


If we write down explicitely H with separate space and time components, we obtain 
1 G, B, — G ; 


1 G, xr 3 
tg E CT COLE 


2mo 


1 
هر‎ gai Et +b,b, -2)= me 
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and finally if we replace G, by i#/c, we get 
G, B, — X Bali 
E p 73 Palic 


0 


le be 1 ch 


۳ CB, BI i 2 
et = E) z ER) A 


LC 


—.— mic? — mo Babs n + by by — 2) 
0 


It is not necessary to derive # explicitely from that equation as we are only con- 
cerned with its derivatives, which can be easily obtained by differentiating the equa- 
tion itself and considering # as a function of the canonical variables x,, t, مره‎ b ” Bu 

1) If the variables x, and t appear nowhere explicitely, the corresponding 
canonical equations 93/2, = - dG,/dt yield 


dre vr 


Thus, in the case when no external force is acting, the particle momentum and its 
energy are constant in time. 
2) Let us now take the derivative with respect to G}: 


2 (o — BP B, ese Jl: ae RR ja + bb.) + 


Aë 8G 
£ 2 G on 
tz(a- ® fa Gb DCE A Gb u) PD 
1 9H Gi Ba ba b» EET | 
ile = he i Ga b,° b, — Ki + b, b, G, i lip Cr bb 


L G 
۳ Da Gh E Ki tarb) (6 — te Be — Gh hb. 


Now, if we take into account the orthogonality and unitarity relations, we see that 
the last term in bracket and the last term on the right-hand side both vanish. So 


we are left with 


0X _ سر‎ Br: Ga Bil — b Cabs 
8G, مه‎ — Ba Ga fih — ba: Caba 


If we remember that 
G, "e Pu N Ga Palha et D G, b, = Mo üp» 


we see that 
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v, being the nonrelativistic velocity dx,/dt. We thus see that we just find again the 
canonical equation 


IK dx? 
8G, ۰ 
3) Let us take now the derivative with respect to b,, the factor mg p, b; being 
dropped, as it vanishes in virtue of the orthogonality relations. So we get 


G 
Ga bı le - SP a. br by) + Gs. (e. - ie, — Cuba b) + 


1 8 ee 1 7 | 
ما‎ a Grey 0 , c ob, 
A By ES E b, b d rs d — my Ba b, 1 b, 


Then, if the orthogonality relations are taken into account, we get 


estu Za Gy by (G, aud Ga Balh§ — bu: G, bj) — m, Brey: b 
9b, Ca — Ba Ga Balh — ba ` Ga Ba 
A very similar calculation provides the derivative with respect to ß,: 
IK _ _ (Go Bulb) (Gu — Bu: Ca Balh — bu: Ga ba) — mo Bu: Ba bal 
Ipu G, — B4: Ga Balh — ba Ga ba 


It will be remembered that in these expressions C, = i #/c is constant in time. As 


for the factors f; b, = — B 1 03, which contain proper time derivatives, they can be 
easily expressed in terms of db,/dt and df,/dt by introducing æ = (1 — v?[c?) * if 
one remembers that 


Ba = فک‎ - © (6, — Ba: 6۸/۱۸8 — ba: Cub): PA. 


This provides finally 


on 3 E BI BR br UA 
LAU CN ES PP OE E 

Han EIU o E MN MN in 
EE Part = 
SB, H 6-0۰6: یط‎ de 


By replacing G, — و60۰6‎ — G,b,:b, by mo ۵ v, the Hamilton equations 


become 
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qe yer Dy pu de 


ee Ber s 
dbp Diet 
=G, ee Ed (LB — گت‎ 
P he de P hà d 
That is 
dë 
T hb) = Gabi vu, 
db; à Ba Ba Ga Ba 
Gr Ve Ae - = Ra: 


which gives just again, after multiplication by a, the covariant equations (5) which 
are equivalent to the second Weyssenhoff equation. 


IV. Conclusion 


In the present paper we have limited ourselves to very simple cases, as well 
with regard to the external field, which was assumed to be zero, as concerning the 
form of the Lagrangian, which was chosen to be linear. Nevertheless our formalism 
can be very useful for several purposes. For instance, we can attempt to frame a fluid 
endowed with intrinsic rotation with an underlying particle structure admitting of 
a chaotic motion. We can approach this problem with the methods of statistical mecha- 
nics which are known to require a Hamiltonian formalism. The above non-covariant 
Hamiltonian will be convenient in the framework of classical statistics, based on the 
use of a definite frame of reference. But it would be perhaps possible to remain in 
a covariant formalism by using a generalized phase space with the variables x,, G,, 
b, , Pu: For this purpose our covariant Hamiltonian formalism may prove very useful. 

On the other hand, by passing from the present classical formalism to a quantal 
one by means of the usual quantization methods, we can put into evidence söme 
new wave equations connected with the spin concept. 


I am indebted to Mr. J. P. Vigier and to Mr. T. Takabayasi for many helpful 


discussions. 
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LETTERS TO THE EDITOR 


IN—PILE HALL COEFFICIENT MEASUREMENTS OF GERMA- 
NIUM BOMBARDED BY FAST NEUTRONS 


Bv B. Buras* J. SUWALSKI 


( Received Jannuary 18, 1960) 


Institute of Nuclear Research, Polish Academy of Sciences, Warsaw 


Lattice defects produced by fast neutron bombardment introduce localized 
energy states into the forbidden energy band of germanium. This causes the change of 
n-type germanium into p-type germanium when irradiated by fast neutrons. In order 
to investigate the nature of this fast neutron bombardment in the neighbourhood 
of the intrinsic region, a semi-continuous in-pile Hall coefficient measurement was 
made. A small electromagnet (placed in an experimental hole of our swimming-pool 
type reactor) and an automatic recording measurement facility were used, and both 
the Hall coefficient and conductivity were measured. The results obtained can be 
seen in Fig. 1. The initial conductivity of the n-type samples was about 0.5 ohm! cm’, 
and the final one of the p-type obtained after 20 hours of irradiation (in a fast neutron 
flux of about 3.1019 cm-? sec!) was about 0.15 ohm-! cm7. Because a small dose of 
fast neutrons was used, it was reasonable to assume that no measurable changes in 
mobilities and scattering factors were caused. On this assumption, the Hall coefficient 
curve was calculated on the basis of a two-carrier model from the conductivity curve, 
provided that the ratio of the electron to hole mobility was 2.05. The obtained discre- 
pancy can be seen in Fig. 2. 

A fairly good fit can be obtained by introducing a three-carrier model. The pre- 
sence in germanium of two types of holes of different effective masses has been indi- 
cated by cyclotron resonance experiments! and verified by magnetoresistance and Hall 
measurements in p-type germanium.? Rencently Walton and Moos? indicated that the 
assumption of 2%, of light holes leads to a mobility of these holes which is six times 
greater than that of ordinary ones. In the described experiment the good fit shown 


ee Tn) re N 


* At present with the International Atomic Energy Agency, Division of Research, Vienna I, Austria. 
(115) 
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in Fig. 2 was obtained with the same percentage of light holes, and the light to ordinary 
holes mobility ratio equals 6.5. In this way we get by a different experiment an addi- 
tional argument in favour of the results obtained by the authors mentioned above 
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Fig. 1. Electrical conductivity and Hall coeffeient R vs. time of fast neutron irradiation 


(see Table I). Some measurements of the Hall coefficient in different magnetic fields 
also support these assumptions. Further work is under way. 


The authors are deeply indebted to Dr H. Rzewuski of this laboratory for discus- 
sion and oriticism. 


TABLE 1 
nee T Percentage Light to ordinary 
go of light holes holes mobility ratio 
Dresselhaus, Kipp, and Kittel! 4? K = T5 
Willardson, Harman, and Beer? 300°K 2% 8 
Walton and Moos? 300°K 2% 


Ours 300°K 2% 6.5 
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Fig. 2. Hall coeffieient vs. time of fast neutron irradiation: dots — experimental points, squares — cal- 
culated on the basis of a two — carrier model, circles — calculated on the basis of a three — carrier 
model 
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STATISTICAL SIGNIFICANCE OF DOUBLE MAXIMUM ANGULAR 
DISTRIBUTION IN HIGH ENERGY JETS 


By J. GreruLA, M. Miesowicz AND P. ZIELIŃSKI 


Institute of Nuclear Research Kraków, and Warszawa 
General Physics Departament, School of Mining and Metallurgy, Kraków 


( Received February 22, 1960) 


In previous papers of our group! è 3 it was found that the angular distributions 
of secondary particles emitted in collisions of high energy nucleons (> 10!? eV) 
with nucleons or nuclei showed the characteristic shape: in coordinates dN/dx vs xi, 
where x; = log tan 0; (©; are the angles of the secondary particles with the primary 
direction) the distribution have two symmetrical maxima. On the contrary the hydro- 
dynamical theories predicted Gaussian distributions, both for nucleon-nucleon and 
nucleon-nucleus collisions. In order to explain the observed shape of the distributions 
we introduced the two-centre model, as we have called it. (See also*,5). In recent 
times several authors have discussed this model®. Since the introduction of the two- 
-centre model was mainly based on the existence of double maximum distributions and 
since on the other hand there were some objections concerning the statistical signi- 
ficance of this fact", we decided to report here the results of a statistical analysis of 
the angular distributions. This analysis was performed on the full available collection 
of jets generated in photographic emulsions by single nucleons of energy higher 
than 101? eV. The total number of jets was 56*. 

The hydrodynamical theories for both nucleon-nucleon and nucleon-nucleus 
collisions predict an increase of the anisotropy (in terms of the ‘dispersion o) with 
increasing energy, the shape of the distribution remaining a Gaussian one, indepen- 
dently of the energy. Therefore we could normalize all distributions to one Gaussian 
distribution. We performed this normalization for all 56 events and the results are 
shown in Fig. 1 as a differential histogram. The intervals on the x-axis correspond 
to equal areas for the Gaussian distribution. For this distribution we obtain a straight 
horizontal line. The deviation from the Gaussian shape is clearly visible on the figure 


* Jets from our laboratories, from published papers, and from different laboratories. We are very 


much obliged to Professor N. Dobrotin for kindly sending us the angular distributions of jets from 
collaborating laboratories and to Dr. D. H. Perkins for the material from H. H. Wills Laboratory, Bristol. 
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